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an 85 m tower, whose elevator carried the traveling microphone. The general 
conclusions drawn from the experimental results were as follows. The 
Inhomogeneities caused significant fluctuations in the instantaneous acoustic 
signal, tut with sufficient time averaging of the measured harmonic levels, 
the results we.’e comparable to resu '■.s expected for propagation in a quiet 
m ^ium.*) (U) 

^Propagation data for the fundamental of the siren approached within 1 dB of 
the weak shock saturation levels. Extra atteniiation on the order of 8 dB 
was observed. The measurements generally confirmed the predictions of 
several theoretical models. The maximum propagation distance was 36 m. The 
narrowbeam arrays were much weaker sources. Nonlinear propagation distortion 
was produced, but the isaximum value of extra attenuation measured was 1.3 dB. 
The maximum propagation distance was J6 m. The behavior of the asymmetric 
vavefoms received in one experiment qualitatively suggested that beam type 
diffraction effects were present.^ (u) 
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FOBEWORD 


This report is an adaptation of Mark A. Theobald's thesis 
"Experimental Study of Outdoor Propagation of Spherically Spreading 
Periodic Acoustic Waves of Finite Amplitude," which was written for the 
degree of Master of Science in Engineering at The University of Texas 
at Austin. Mr. Theobald was enrolled in the Mechanical Engineering 
Department, and his degree was granted in May 1977* The work was 
carried out at Applied Research Laboratories. 

The research, which began in 1975> represents the first phase in 
a study of the propagation of finite amplitude sound outdoors. In this 
phase, periodic sound is used. Random noise is to be used in the second 
phase. Jointly carrying out the first phase with Mr. Theobald was 
Don A. Webster. Some of the results not described in this report may 
be found in the appendices of Mr. Webster's report, "Saturation of Plane 
Acoustic Waves and Notes on the Propagation of Finite -Amplitude Spherical 
Waves," Applied Research Laboratories Technical Report ARL-TR-77-^ 
(January 1977)* 

The research was supported by the National Aeronautics and Space 
Administration under Contract NASl-l^l60, by the Air Force Office of 
Scientific Research under Contract F^^20-76-C-004o, by the Office of 
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Oceanic and Atmospheric Administration under Grant 0k-^~022-12. 
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ABSTRACT 


The outdoor propagation of spherically spreading sound waves 
of finite amplitude was investigated. The main purpose of 'he experi- 
ments was to determine the extent to which the outdoor environment, 
mainly random inhomogeneity of the medium, affects finite amplitude 
propagation. Periodic sources with fundamental frequencies in the range 
6 to 8 kHz and source levels SPL, from 1^ to l4g dB were used. The 
sources were an array of 7 to 10 horn drivers and a siren. The propaga- 
tion path was vertical and p6'.rallel to an 8*5 m tower, whose elevator 
carried the traveling microphone. Measurements were made in late summer 
under a wide range of weather conditions. The ground level wind speeds 
were from 0 to 24 km/h. Measurement error was of the order 0,5 dB for 
calm air. 

Several theoretical approaches were used to explain the data. 

A graphical method is used to gauge the expected importance of nonlinear 
effects, k computer algorithm is described for use with finite amplitude 
spherical waves in lossy media. Several ad hoc propagation models are 
applied. 

The general conclusions drawn from the experimental results 
were as follows. The inhomogeneities caused significant fluctuations 
In the instantaneous acoustic signal, but with sufficiem; time averaging 
of the measured harmonic levels, the results were comparable to results 
expected for propagation in a quiet medium. 


Vi 


Propagation data for the fundamental of the siren approached 
within 1 dB of the weak shock saturation levels. Extra attenuation on 


vii 


the order of 8 dB was measured. The measurements generally confirmed the 
predictions of several theoretical models. The* maximum propagation 
distance was 56 tn. A power output at the fundamental frequency of kio W 
was f*stimated. 

Experiments using the narrowheam arrays showed that nonlinear 
propagation distortion was produced. The maximum value of extra attenua- 
tion measured was 1.5 dB. This value was approximately 1-5 'iS less than 
predicted. The maximum propagation distance was 76 m. The hehavior of 
the asymmetric waveforms received in one experiment qualitatively 
suggested that beam type diffraction effects were present. 

The role of diffraction of high intensity sound wa\as in 
radiation from a single horn was briefly investigated. 
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CHAPTER I 


INTRODUCTION 

The general research topic of this report is the propagation of 
finite amplitude airborne sovind in the outdoor environment. Jet and 
rocket engines, explosions, and atmospheric soimders ("acoustic radar") 
are common sovirces that may produce sound intense enou^ to cause finite 
amplitude propagation distortion to he important. Only a few studies of 
spherical periodic sound waves of finite amplitude in air have been 
ceniTied out in the laboratory; even fewer have been done outdoors. 
(Underwater measurements have been made in the fiel:^ as well as in the 
laboratory.) Thus, because oxir ignorance is great, the need for outdoor 
studies is great. Two sample problems miy be mentioned, (l) The attenua- 
tion of broadband noise from jet aircraft cannot now be precisely 
calculated. The effective attenuation of the hi£^ frequency spectral com- 
ponents of jet noise is considerably less than expected when only spheri- 
cal sjjreading and atmospheric absorption are accounted for. The lower 

1 ** 

attenuation at hlg^ frequencies might be due to nonlinear effects. 

( 2 ) Nonlinear losses may limit the range and effectiveness of atmospheric 
sounding Just as they do in underwater sounding. On the other hand, non- 
linear effects may not be as important outdoors as they are in the 

Thf atmospheric acoustic sounder has been operational since I968 (see, for 
example. Ref. 5). The device is used to detect atmospheric inhomogene- 
ities, such as wind sheen: layers and thermal inversions, by receiving sound 
pulses reflected from the inhomogeneities. The acoustic source is a group 
of compression horn drivers feeding a common horn or parabolic reflector; 
an alternative source is the arrangement of the drivers in a plane array. 

References are listed in ntimerlcal order following the appendices. 
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laboratory. It has been speculated, for example, that inhomogeneity of 
the outdoor medium might neutralize or inhibit nonlinear distortion 
effects. In any case, very high sound levels are encountered ever more 
frequently. It is time to replace speculation about the role of non- 
linearity vith analysis and finn experimental evidence. 

The reason for performing outdoor rather than laboratory 
measurements is, of course, that the medium is different. Turbulence, 
wind, and thermal gradients are not normally encountered in the labora- 
tory; we do not know very much about their effects on the propeigation of 
finite amplitude waves. Our experiments are sm initial step toward 
experimentation in an uncontrolled environment. 

Our objectives were to measure the finite amplitude distortion 
of a high intensity audio frequency sound wave outdoors, to compare the 
measurements with existing theoretical predictions, and to record any 
oth . nonlinear phenomena encountered. To simplify the interpretation of 
our results, we kept the experiment as simple as possible; thus, a verti- 
cal propagation path was chosen to avoid complications from ground 
reflections- The experiments were restricted to conditions of low wind 
speeds to keep random medium effects to a minimum. Two different sound 
sources were used: an array of horn drivers and a siren. 

Our research was motivated partly by two very practical 
questions, (l) Are nonlinear effects important in the propagation of Jet 
noise? Because the experiment needed to be kept simple, we decided to 
use periodic waves, not noise. Since the results obtained in the present 
study were positive, that is, that nonlinear effects can be important in 
outdoor propagation, the next step will be to employ a source that emits 
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intense noise. (2) Are nonlinear effects importemt for acoustic sounding 
of the atmosphere? Because arrays of horn drivers simileir to ours have 
been used for acoustic sounding in this application, our results may 
indicate a partial answer to this question. 

The report is divided into the following chapters. 


I. 

Introduction 

II. 

Review of the Literature 

III. 

Application of the Theory 

IV. 

Experiment Design 

V. 

Array Experiments 

VI. 

Siren Experiments 

VII. 

Conclusions 


Appendices 


The reader familiar with theory concerning the propagation of finite 
amplitude spherical waves may omit chapter II. The considerations in 
directly applying the theoretical models to our experiments and a com- 
puter algorithm applicable to spherical wave propagation are described in 
chapter III. General criteria for the design of o\or propagation experi- 
ment are given in chapter IV. The evolution of the acoustic sovirces is 
described in chapters V and VI. The experimental methods are explained 
and our results are compared with the existing theory. Chapter VII 
contains closing remarks. Finally, the appendices contain material of a 
specialized natiare. The performance and construction details of an 
anechulc chamber for use at hlg^ audio frequencies are given in Appendix A. 
The role of diffraction in hl^ amplitude radiation from a horn is 
discussed in Appendix B. 



CHAPTER II 


REVIEW OF THE LITERATURE 


A. Introduction 

In this chapter we consider existing theoretical and experimental 
results relevant to the research problem. We seek to describe the distor- 
tion and extra attenuation suffered by a directional spherical wave 
propagating in an outdoor environment. The extra attenuation is the 
diminution of the fundamental over and above which a small signal wave 
would suffer. In this report, the terms "small signal" and "linear 
theory" denote both spherical spreading and atmospheric absorption. Any 
other atmospheric effects are mentioned separately. There are several 
analytical as well as ad hoc models available to describe nonlinear propa- 
gation under various conditions, but we shall limit our discussion to the 
models that were applied to the experiments. 

In most theoretical treatments of the propagation of finite 
amplitude waves, a perfectly sinusoidal source wave is assvimed. In 
practice, source distortion or effective source distortion compromises 
this assumption. It is therefore necessary to find ways of modifying or 
applying the theoretical results so that more realistic source waveforms 
may be considered. This question, however, is postponed until the next 
Chapter. For now only' purely sinusoidal source waveforms are considered. 

B. Theoretical Literature 

1. Propagation in a Homogeneous, Quiet Medium 

We begin the review with a qualitative description of the 


distortion of an initially sinusoidal wave propagating in a homogeneous, 

4 
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quiet medium (see, for example. Ref. 4). We assume that ordinary 
absorption is negligible. To fixrther simplify the discussion, we 
initially confine ourselves to plane wave propagation. First, it is 
helpful to define the term "wavelet." A wavelet is a point on a waveform} 
a specific value of particle velocity u or acoustic pressure p is asso- 
ciated with this point. Each wavelet propagates with a speed depending 
on its particle velocity, as expressed in the foarmula 

= c + Pu . (2-1) 

u=constant 

Here, x is the position, t is time, c^ is the small signal sound speed, 
and p [=(7+1)/2 for air, where v is the ratio of specific heats] is the 
coefficient of nonlinearity. The physical mechanisms that lead to the 
second term in Eq. 2-1 are convection and the nonlinearity of the medium’s 
pressure density relation. For a small signal wave (u«c^), all wavelets 
travel at the same speed. But for increasing values of u, the propagation 
speed of the wavelets departs from consteuicy. In particular we see from 
Eq. 2-1 that the wave peaks (positive values of u) travel faster than the 
troughs (negative values of u) . Oxir sinusoidal wave will therefore 
distort as it propagates. At some distance x, a shock or discontinuity 
will form in the wave. The unrestricted application of Eq. 2-1 beyond 
this distance leads to difficulty because a multivalued wavefoim is pre- 
dicted. Because such a waveform is physically disallowed for acoustic 
waves, Eq. 2-1 must be modified. Before taking up the modification, 
however, let us first examine the wave equation associated with Eq. 2-1. 


dx 

dt 
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fluid ir. 


The epp'^'OxiiQate wave equation for pla'»t 'ives in a lossless 


u 


X 



uu 


t ‘ 


( 2 - 2 ) 


where t' is the time t-x/'i . ...ict equation is less 

C 

convenient for our purposes and regli, . more accurate, even up to sound 
pressure levels of the order 165 dB) . "'he solution of this equation 
satisfying the boundary (source) condition 


u(0,t) = g(t) , 


(where g(t)=u^ sin cut for the present discussion) is 


u = g(<p) 


Cp = t' 



(2-5) 


The notation is as follows: qp is the time at which a particular wavelet 

leaves the source (i.e., t=cp at x=0) and u^ is the peak particle velocity 
at the source. Equation 2-5 is called the Earnshaw solution. (Actually, 
it is an approximation of the solution given by Earnshaw.^) It is valid 
for distances x<x, where x is the shock formation distance. A sinusoidal 
wave forms a shock at x=l/Bek, *fhere e[=u^/c^] is the acoustic Mach number 
at the source and k is the wave number. At this point it is convenient to 
introduce the distortion distance variable cr, which is defined for plane 
waves as x/x. Thus for a wave from a sinusoidal source, o=Bekx, and a=l 
signifies shock formation. 


Throughout this work, sound pressure levels are referred to 2xl0"^ N/m^. 
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We now consider the question of modifying Eqs. 2-1 a>'id <^-5 so 
that wave trains containing shocks may be described, fhe Earnsnaw 
solution (Eq. 2-3) may still be applied to the continuous portions of the 
waveform between shocks. The Rankine-Hugoniot shock relations, vi uich 
express conservation laws for shock propagation, may be combined to jield 
the useful relation 


dx 

dt 


shock 


- =0 ♦ I 


(2-4) 


where u^ and u^^ are the particle velocities immediately ahead of and 
behind the shock, respectively. The expression indicates that the velocity 
of the shock is approximately the mean of the .ravelet speeds immediately 
ahead of and behind the shock. Equation 2 -^ then replaces Eq. 2-1 in 
describing the propagation of the shock front. Descriptions of the 
merging of shock waves are provided by Refs. 4 and 5* The •omblnation of 
the Earnshaw solution with Eq. 2-4 is termed "weak shock theory." Few 
Initial waveforms are simple enough for analytical sol.utions of the weak 
shock theory equations to be found. On the other hand, the equations 
themselves are conceptually simple and amenable to solution by computer 
methods. Pestorius^ has developed a computer aigorlthra that may be used 
for an arbitrary sorxrce waveform 

In the propagation of an initially sinusoidal wave, the wave 
will distort to form a full-fledged sawtooth at the position x=3x(o=3).^ 

The sawtooth is the most nearly stable periodic waveform for a finite 
eunplltude wave. Dissipation occurs at the shocks in the sawtooth. 
Simultaneously, energy from the fundamental is pumped into the higher 
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harmonics. This process keeps the wave a sawtooth. The sawtooth 
waveshape is approximately maintained out to a distance x =l/a-x, 

(^ thay )^ where a is the attenuation coefficient for ordinary absorption in 
the medium. From approximately this point on, the wave becomes a quasi- 
slnvisoid. Linear theory is then roughly sufficient to predict the 
propagation. 

The major weakness encountered with weak shock theory in our 
application is that the only dissipation included is caused by irreversi- 
ble processes at the shocks. Ordinaiy small signal absorption is ignored. 
In practice, there are many instances in which ordinary absorption is not 
negligible for waves containing weak shocks. In fact, absorption may 
prevent shock and/or sawtooth formation by damping the wave very quickly. 
Weak shock theory should not be applied to waves too weak to form shocks 
(see, for example. Ref. ^). Tlius the problem of simultaneously accoiuiting 
for both nonlinear and ordinary absorptive losses must be considered and 
will be discussed in more detail ]ater in the chapter. 

The application of weak shock theory is restricted to pressure 
levels for which the underlying appreximations are valid. The upper limit 

n 

for plane waves in air is a sound pressure level of approximately I 65 dB. 
This restriction was of minor importance in our experiments. 

Turning now to the propagation of diverging spherical waves, we 
note that the approximate (lossless) wave equation analogous to Eq. 2-2 is 



where the retarded time Is now t'=t-(r-r )/c and r is the sovirce radius. 

' o" o o 

The nonlinear distortion process slows down dramatically for spherical 
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waves, as compared to the distortion process for plane waves. The reason 
for the slackening is the reduction in amplitude caused by spherical 
spreading. The mechanisms causing the distortion are the same as the 
mechanisms for plane waves, but their effect is lessened by spherical 
sx«readlng. The Eamshaw solution, Eq. 2-5, and weak shock theory may be 
applied to spherical propagation by making the replacements u-»ru/r^ and 
x->r^ln(r/r^) . For spherical waves from a sinusoidal source, the 
expression for o is 

0 = pekr^ In^f-^ , 

which shows that a increases much more slowly with distance thcui it does 
for pleme waves. The inverse relation is 


r 


a/pskr 

e 


The values o=l and c=5 still Indicate shock and sawtooth formation, 
respectively. For example, shocks form at the distance 

l/p€kr^ x/r^ 

r = r e = r e 

o 0 

Therefore, if the equivalent plane wave shock formation distance x is much 
greater than the source radius, the spherical wave must travel a very 
great dlsteuice to distort into a shock wave. As in the case of plane 
waves, there is also a distance r (a ) beyond which nonlinear effects 
become approximately negligible. The formula for r^ ^ is given in 
chapter III. Directional sources may be treated for the farfield case by 
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a procedure described in Ref. 9* jecause our measurements were made along 
the axis of the sources, we ignore directionality in the following 
material. 

Included in the approximations leading to Eq. 2-9 is the 
farfield assumption, that is, kr»l. In the farfield of a spherical 
source, the peurticle velocity and the acoustic pressure p are nearly in 
phase. Althou^ the solution of Eq. 2-5 is most natinrally written in 
terms of u, the condenser microphones used in our exx>erifflent8 measure 
pressure. The small signal Impedance relation 

P = PqC^u (2-6) 

may be used to relate p and u for the sound levels amenable to solution 
by weak shock theory.^ 

In the preceding discussion, nonlineeu: distortion was viewed in 
the time domain. The view from the frequency domain is often equally 
important. In other words, given a source waveform, what is the spectrum 
of the wave at a certain propagation distance? We consider a solution to 
this question in terms of the Fourier series 

P = — ^ sin n[ojb - k(r-r^)] , (2-7) 

subject to the boundary condition 

p = Pj^Q sin 03b at r=r^ 

We shall discuss several approaches for calculating the harmonic ampli- 
tudes B^. We first define the often used abbreviation EXDD. The extra 
n 
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attenuation of the fundamental d.ue to nonlinear effects- may be expressed 
as the ratio (in decibels) of the fundamental amplitude to the combin^'d 
effects of spherical spreading and small signal absorption. This quantity 
is denoted by EXDB, where the acronym stands for the "extra attenuation 
expressed in decibels." 

Perhaps the simplest expressions for the harmonic amplitudes are 
the expressions for the amplitudes of strong waves, that is, waves in 
which ordineury absorption compared to shock dissipation is negligible. 

The Pubini solution (see, for example. Ref. is valid in the preshock 
region (c^l). In this case, the harmonic amplitudes are given by 

B = -§- J (na) , (2-8) 

n ncr n' ' ' ' 

where J is the nth order Bessel function of the first kind. In the 
n 

sawtooth region (3^og g^^^ ) , the weak shock result is 

\ ■ ktIst • 

(We note for later reference that the second harmonic of a sawtooth wave 
is 6 dB below the fundamental.) The expression valid in the region l<o<5 
is more complicated euid is given i:. Ref. 8. 

The analysis of moderate and weak waves musu Include the effects 
of ordinary absorption. Several models have been proposed to simultane- 
ously consider nonlinear distortion and absorption. Burgers' equation 
describes simple wave flow in a thermovlscous fluid (see, for example. 

Kef. 4). No exact solutions are known for spherical waves, but some 
approximate solutions have been fomd. For example, a perturbation 


I 
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solution has been obtained for the range prior to shock fon.ation.^*^^^ 

In describing the results of our first array experiments, we shall make 

use of the second order perturbation solution. Natigol'nykh et al. found 

an approximate solution of Burgers* equation for the sawtooth wave 
12 15 

region. Tepper approximated Na\igol'nykh*s solution by a Fourier 
13 

series, and Cary used this solution to calculate the extra attenuation 
(EXHB) of the fundamental. In chapter VI we compare values predicted bv 
Oary with values of the extra attenuation measiired during our siren 


experiments. 

6 l4 

Pestorius and Anderson have programmed weak shock theory for 
the propagation of plane and spherical waves, respectively. Pestorius 
corrected the peurticle velocity values for tube wall attenuation and dis- 
persion during the distortion process. The calculations were in good 
agreement with experiments performed in a plane wave tube. Anderson 
generated computer solutions for the propagation of a spherical N wave in 
air. The solution was improved by including air absorption in the calcu- 
lations. Anderson's computations closely matched his experimental 
results. To predict received waveforms for our experiments, we modified 
Pestorius* program by including spherical spreading and atmospheric 
absorption. 

In addition to the models given above, several ad hoc models 
have been proposed to simultemeously account for both absorption and non- 
linearity. These models are ad hoc in the sense that they sure not direct 
attempts at solving the wave equation. Safar'*'^ and Pemet and Payne 
present models valid in the preshock region for absorption of arbitrary 


5 
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ftrequency dependence. The sound generated by our sources was generally 
too Intense to apply these two models. 

Merklinger et al.^*^ used a method suggested by Westervelt^^ to 
calculate the extra attenuation for spherical waves. Westervelt ass\jmed 
that energy from the fundamental is used in second harmonic generation and 
Hmnli signal absorption. Interactions involving hl^er order heumionlcs 
are Ignored. The suggested model is 


div I = - 2al + 

where I is the intensity of the fundamental^ I=|I|, the angle brackets 
denote a time aerage, smd q is the nonlinear acoustic source density per 
unit volume 


q = 




P c 

o 


(p(« . . ...)‘ 


Here p^^^ is the small signal solution for the fundamental at the retarded 

time t' and p^^^ is a second order correction involving the second har- 

15 

monlc. Using a perturbation solution due to Safar and solving a 
Bernoulli equation for I, Merklinger found 


EXDB = 10 logj^Q 



( 2 - 10 ) 


where a =Sekr euid H is a double exponential integral that •mist be 

O 0 

evalxiated nuznerically. We compare our experimental results with predic- 
tlons based on Eq. 2-10 in chapters V and VI. 
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Rudnick, ^ while Investigating the attenuation of plane sawtooth 
waves in a duct, made the ad hoc assvonption that the sawtooth decay rate 
is simply the sum of the decay rates due to small signal attenuation (tube 
wall effects in his case) emd nonlinear effects. The assumption is 
expressed as 

= - (“tPpea. * 

where is the acoustic ja-essure at the sawtooth peak and is an 

absorption coefficient for the tube wall attenuation. More recently, 
Webster^^ used the Rudnick model to predict the saturation level of the 
fundamental for plane waves in a tube. Webster also presented a solution 
for the fundamental amplitude for the case of spherical waves. In 
chapter VI we compare data from our siren experiments with predic- 
tions based on the Rudnick model. 

2. Effects of Reindom and Stratified Media 

Two basic types of inhomogeneity of the mediiion must be considered 
in planning an outdoor propagation experiment. First, the atmosphere is 
stratified or layered with respect to the fluid parameters. The stratifi- 
cation is especially evident in vertical propagation at heists that are 
nonnegllgible when they are compared to the scale height of the 
atmosihere. (Because we consider only vertical propagation, the efiec s 

* 

The scale hel^t 6 is the basic length scale related to stratification of 
the atmosphere. It is Cf^/{7g) (about 8 km), where g is the acceleration 
due to gravityj G Ir the hei^t of the atmosphere if the density were 
constant. 
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of refraction are neglected here.) Second, rsmdom fluctuations of the 
medium paramet«^rs often accompany wind turbulence and other disturbances. 
While the propagation distance for the experiments described in this 
thesis was very short, future experiments at longer ranges are envisioned. 
It is therefore appropriate to mention the effects of inhomogeneities 
here. 

A sizeabL. body of theoretical literature describing firite 

amplitude propagation in stratified media has evolved. The atmosphere is 

often modeled as an Isothermal medium subject to the force of gravity. 

The ambient density for an isothermal atmosphere decreases exponentially 

with height, but the small signal sound speed remains constant. For an 

analysis of the effects of stratification for this case (and a list of 

20 

other Russian works on the subject), see Romanova. The assumption of 

constant (ambient) temperatiare is not required, however. Carlton and 
21 

Blackstock treated the vertical proxjagatlon of plane waves in a mediian 

where the small signal sound spaed, ambient density, and fluid nonlinearity 

22 

vary with height. Nayfeh described quasiplane wave propagation in 
ducts with varying cross section and fluid properties. Both references 
may be applied to plane propagation in the atmosphere when the functional 
dependence of the fluid parameters is known. The effect of stratification 
is to amplify or reduce the amplitude u^ of the particle velocity relative 
to the sound speed c^. Because finite amplitude distortion depends on 
Uo/Co# stratification has the effect of either speeding up or slowing down 

t 

the distortion. For example, for a wave propagating vertically in the 
atmosphere, the distortion generally accelerates if the wave is traveling 
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upward but decelerates if the wave is traveling downward. The details 
depend ultimately on the specific medium conditions and the soiirce 
amplitude . 

The theory of linear propagation through a meditan containing 
random inhomogeneities is well known (see, for example. Refs. 23 and 24). 
The effect of a random medium on nonlinear propagation is the subject of 
several recent publications. A difficulty in applying random medium 
theory to an experimental situation is that the distribution function in 
time and space for the inhomogenelties is, in general, unknown and very 

26 

difficult to measure. However, Fridman and Pelinovsky state that, for 
finite amplitude propagation, the amplitudes of all harmonics above the 
fundamental are larger thah the amplitudes in the absence of fluctuations. 
The Increase depends on the distribution functions of the parameters. 
Lacking the instrumentation reqvdred to estimate the distribution func- 
tions, we assumed that, by confining our Initial experiments to conditions 

of low wlndspeeds and by averaging the data over a sufficiently long time 
27 

span, we could ignore the effects of random fluctuations on nonlinear 
distortion. We found this assumption to hold true in our experiments. 

5 . Diffraction Effects 

Some recent theoretical work from the USSR concerning the 
effects of diffraction on strong acoustic beams is of interest. The 
theoretical models describe the propagation of a well-collimated circular 
beam of plane waves. The amplitude distribution across the beam is 
assumed to be 

P = P^o 


( 2 - 12 ) 
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where r is the coordinate transverse to the beam and b is the effective 

beam radius. (This assumption is commonly made in raden: theory.) 

28 

Rudenko et al. gauge the relative importance of diffraction to nonlinear 
distortion by the parameter 

N = ^ . (2-13) 

pe(kb)‘" 

As N-»0, nonlinear effects dominate the behavior vhile> for N-»», diffrac- 
tion is most important. Analysis for the lossless case where diffraction 

yo * 

is nonnegligible shows that the beam diverges after x is surpassed. The 
most intense section of the beam (on and near the axis) is damped by 
dissipation at the shocks. The wave becomes asymmetrical^ that is, the 
peaks are sheirpened and the trouts sure rounded. Figure 2-1 shows the on- 
axis solutions for the cases N=0 and N=0.4. Bakhvalov et al. show theo- 
retical waveform variations on and off the axis as functions of soxirce 
amplitude, the distance off-axis, and propagation distance. No solutions 
of the governing equations for arbitrary source waveforms are available. 

The well-collimated beam described above is not encountered in 
practice, but the model may apply to highly directive radiators at large 
distances where the waves are nearly planenr. The received waveforms in 
o\ir array experiments may show qualitative evidence for the "peaked" 
waveform described above. In our discussion of the data we estimate the 
value of N for our experiments. For a discussion of the diffraction of 
hi^ intensity waves from a horn, see Appendix B. 

* 

The definition of the discontinuity buildup length in Ref. 28 corre- 
sponds to the "well-formed sawtooth distance" x(a=3) in our notation. 
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AXIAL WAVEFORMS IN A NARROW BEAM OF FINITE AMPLITUDE 
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C. Ex'nerimental Studies 

There have heen very fev experimental studies of periodic 

* 

spherical waves of finite amplitude in air. The only known outdoor 

54 

measurements were made by Widener and Muir,''^ who investigated the per- 
formance of a parametric array operating in air. (The term "parametric 
an*ay" describes the nonlinear interaction of two collinear single fre- 
quency sound beams of finite amplitude. See, for example. Ref. 55.) The 
difference frequency generated by the interaction was investigated along 
a 20 m long horizontal path over water. While other outdoor experiments 
using very powerful sources have been reported, no attempt was made in 
these experiments to quantify the role of nonlinear effects on the data 
(see, for example. Refs. 36 and 37) • 

The laboratory (indoor) studies most nearly related to our work 

38 39 4o 

are those of Allen and Rudnick, Allen, and Shin. Allen and Rudnick 
measured the beam patterns for a siren generating an acoustic output of 
84 to 176 W and pointed out several interesting chetracteri sties of the 
intense sound field. 

4l 

Allen Investigated the radiation from a St. Clair generator 
operating at l4.6 kHz in a 2 m long anechoic space. The St. Clair 
generator functioned as a plane baffled piston and, as a resonant device, 
generated a nearly pure sinusoidal wave at a maximum SPL of 161 dB. Allen 
recorded such finite amplitude effects (now commonly known) as acoustic 
saturation, blunted beam patterns, and poor side lobe suppression. 
Saturation is the limiting process whereby the received level at a field 

\e duly note the large literature on sonic booms (see, for example. 

Refs. 31 and 32) and explosion sounds (see, for example. Ref. 35). 
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point is independent of the source level. Increases in the source le\'el 
are nearly nullified by extra attenuation. The directivity alterations 
eire caused by the extra attenuation incurred at sections of the sound 
field with hi^er intensity. A major lobe that is normally rounded 
becomes squared off as the amplitude in the center decreases relative to 
the side lobes. 

To study the attenuation of spherical sawtooth waves^ Shin 
utilized a siren and conducted the study in an anechoic chamber. He 
corrected the weak shock decay rpoe for spherical spreading auid obtained 
excellent agreement with his experiment. The data cover the frequency 
range 8.5 to 11 kHz for extrapolated source levels of I59 to I*!-? dL at 
1 m, over distances as great as 1 m. 
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CHAl'TPER III 

APPLICATIOH OP THE THEORY 

A. Effective Source Parameters 

One of the difficulties in applying the model 3 described in 

chapter II is that the source conditions assumed in the models are usually 

different xrom the conditions realized in the experiments. -Except for the 

computer al 7 ,orithm described in section B.l. ^ this chapter, the models 

used here are restricted to sound that is a pure sine wave at a reference 

point r^ in the farfield. In practice, however, the wave is usually 

already distorted when it reaches the farfield. There eire two reasons 

for such distortion, (l) There is inherent distortion in the transducer. 

Inlierent distortion was particularly strong in our siren soui*ce. (2) The 

signal is subject to finite amplitude distortion as it travels through 

the nearfield. Our method of dealing with "source distortion" is to 

measure the distorted waveform at an Initial point r=r^ just outside the 

nearfield and to extrapolate baclarard until we reac.h an effective source 

* 

radius r^ at which the waveform is a pure slnus^ia. Ine dimensionless 

amplitude at this point is denoted e. 

Our estimates of r and e eure based on a measurement of the 
o 

sound pressure levels of the fundamental and second harmonic at the 
initial point r^ . We use the second harmonic rather than other harmonl cs 
because, in our experiments, more energy was present there than In the 

For piston sources operating underwater, it has been found empirically 
that a value of r^ between Rq/J and acceptable. Here R^ is the 

Ra^'lelgh distance given by the piston area divided by the wavelength. 

See, for example. Ref. 2. 
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higjaer harmonjcs. Consider again Eq. 2-7 and let denote the pressure 
amplitucLc the nth ha'^monic. We write 


(4 

f 

i 




^>1 = 


^ o^>10 


B, 


and 


Pg = 


’^o^lO 


B. 


; I 


where B^ is given hy the appropriate formula (for example, Eq. 2-8 or 
2-9), depending on whether shocks are present and at what stage of 
development. A measurement of p^^ and Pg is made at the distance r^. We 
are left, therefore, wit’i a system of two equations in two variables (r^ 
and Pj^q)* In principle, we may solve the system to yield the desired 
source parameters. For example, if r^^ is in the preshock region, that is, 
o^<l, we have 


2_ 

'i ^i 


Pt . = — — J (o.) 

^li Tj a. 1' X 


(5-1) 


V * 

4 ' T 

>> 


and 


_ ^o^IO 1 T ^ 
^2i ■“ r, 0. 


(5-2) 


Forming the ratio of Pg^ to p^^, we obtain 


X 1 

} 


Pll 2JjT^ 


(5-5) 


from which the value of 0 ^ may be found. The smallest value of 0 ^ satis- 
fying Eq. 5-5 ds selected. Equation ^-1 may then be solved for r^p^^Q. 


I 
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Using this value emd the expression for [=p€kr^ln(rj^/r^) ], we obtain 
the source radius 


= 


-Oi/pekTo 


(5-M 


Finally, (and therefore e) is computed from the value 

shocks are present, one must use the sawtooth or transition expressions 

for p^ and Pg in order to find 

In the extrapolation procedure described here, it is tacitly 
assuiLsd that the fundamental and second harmonic are in phase. In 
general, this condition is not guaranteed to hold. In addition, a source 
may produce an asymmetrical waveform. We are therefore justified in the 
extrapolation procedure only to the extent that agreement with the experi- 
ment warrants. 

8 

The range i>arameters for our experiments are calculated by 



' o 

(5-5a) 



(3-5b) 

euid 

/o \ r /r \1"^ 




(5-5c) 


where c^=pekr^. Equation 5-5c is transcendental and must be solved by 
iteration. Range parameters based on an effective source are listed with 
the propagation curves (shewn in chapters V anc VI). Ihe parameters are 
subscripted "eff" to suggest the effective source. 


2h 


The calcxilation of extra attenuation based on the effective 
80\jTce parameters also warrcmts a word of explanation. Measured extra 
attenuation must refer to the level of the fundamental at the source posi- 
tion r^, that is, the point from which the effective source parameters 
were extrapolated. The predicted va?.ue of extra attenuation for our 
source is then the difference between the values predicted by the models 
for the disteuace r(r>r^) and r^. 

B. Theory Applicable to Both Moderate and Strong Waves 
1. Computer Prope^tatlon by Weak Shock Theory 

The propagation of a spherical wave of finite amplitude in a 
lossy medium may be described by means of a computer program. We modify 
Pestorius' weak shock propagation algorithm for plane waves^ by correcting 
the computations for spherical spreading. For a loss mechanism, we sub- 
stitute atmospheric absorption in place of the boundary layer attenuation 
used by Pestorius. D. Watson of Applied Research Laboratories provided 
additional, more flexible input and output schemes. Only the basic 
steps in the computations are outlined below. See Ref. 6 for a detailed 
discussion of the techniques and algorithms used to perform the calcu- 
lations. 

The propagation of a finite amplitude spherical wave consists of 
three simultaneous processes: (l) spherical spreading, (2) nonlinear 

'’stortion, and (3) ordinary absorption. For a sufficiently small propa- 
gation distance, one may consider the processes to be Independent of each 
other and may calculate their effects sequentially (see, for example. 

Ref. 6) . While a process may be described in either the time or frequency 
domain, the functional representation in one domain is often simpler. 
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The cconputer's forward and inverse fast Fourier transforms (FPT and FFT**^, 
respectively) link the two domains. Because of the restriction kr»l, we 
are Justified in using Eq,. 8-6 to relate particle velocity euad pressure. 

The flowchart in Pig. 5-1 provides an overall view of the 
computer algorithm. The functional form for each part of the propagation 
process is listed below. 

(l) The particle velocity decreases by spherical spreading at 
a rate of 6 dB per doubling of distance. This relation ie expressed in 
the time domain as 


= rl £ar 


(5-6) 


where Ar is the size of the propagation step. 

( 2 ) The finite amplitude distortion is viewed in the time 
domain as an amplitude dependent shift in the arrival time for a particu- 
lar wavelet. If the wavelet were of infinitesimal amplitude or were a 
zero of the waveform, the time required to traverse a distance Ar would, be 


t 

° *^0 


But for a finite amplitude wavelet whose particle velocity is u, the 
travel time for the same distance is 


t =t - ^ 

f c + Pu o 2 . Q 

° Co + c^pu 


To first order in e (where u«c^), we may drop the term c^Pu in the above 
equation to obtain 
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(3-7) 


The Euialogous form for the arrival time t of a weak shock front is 

6 

(see Eq. 2-4) 


t 

s 



p(u^+u^)Ar 


2c 


(3-8) 


The computations required for the merging of shocks are given in 
Ref. 6, Our version of the ccmputer program is not included in this 
report because the published version by Pestorius may be easily modified 
by including spherical spreading (Eq. 3-6) and by replacing tube wall 
losses by atmosiiieric absorption. Because nonlinear distortion depends on 
the ratio of the particle velocity to the small signal soxmd speed and our 
measurements were taken at widely varying temperatures, we were careful to 
use the appropriate speed of sound in the calculations. 

( 3 ) l^e wave is damped by atmospheric absorption. The damping 
increases with frequency and the calculation is performed in the frequency 
domain. The functional form for the atmospheric absorption alone is 

-a At 

U^(r+Ar) = Uj^(r) e " , (3-9) 

where is the magnitude of the nth harmonic in the Povirler series for u 
and is the respective absorption coefficient. Dispersion need not be 

An appropriate place to insert the calculation for spherical spreading in 
the published program is following step c4o, p. 173 of Ref. 6. The accumu- 
lated propagation distance prior to calling the subroutine is required for 
the calculation. 
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considered, because in atmospheric propagation it is negligibly small at 
normally encoxintered frequencies. 

In our experiments atmospheric absorption is a much weaker 
effect them spherict.1 spreading and is often weaker than nonlinear dis- 
tortion. The absorption calculation occurs in the frequency domain while 
both the distortion and spreading calculations are made in the time domain 
(see Pig. 5-1) • Dius, computer time is saved (the FPT operation is com- 
paratively slow) hy correcting for absorption at larger distance incre- 
ments than those used for distortion and spreading. The absorption 
coefficients eure calculated from a formula given by Bass. The resulting 
coefficients eure accurate to ±5^ over the temperature range 275 to 515 ‘*K, 
the relative humidity range 0 to 100^, and the frequency range 50 Hz to 

10 MHz at a pressiire of 1 atm. The addition of an algorithm based on 

l4 

Bass' formula is an improvement over the procediire used by Anderson, who 
assumed that the coefficients simply increase with the square of the fre- 
quency. The improvement should be significant at frequencies where 
relaxation effects are important. 

By repeatedly applying Eqs. 5-6, 5-9» and 5-T or 5-8 to every 
wavelet, one may simulate the propagation of the wave. 

The input waveform to the program was taken from a measurement 
at r=r^^. The point r^^ was in the far fie Id of the source but was as close 
to the source as feasible. A Hewlett-Packard 981 OA programmable calcu- 
lator with a s86kA digitizer was used to digitize the oscillogram. The 
minimum resolution was 0.25 mm. The waveform was smoothed by the 
calculator to minimize digitizing Jitter. The resulting data were 
filtered to remove any dc pressure bias (do bias would cause T'roblems 
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dijring program execution) and were scaled to match the time-averaged level 
of the fundamental. Discussion of further details related to program 
execution is delayed until chapter V. 

2. Merkllnger, Mellen> and Moffett Model 

In order to apply the method developed hy Merklinger et al. for 
the calculation of the extra attenuation of the fundamental (Eq, 2-10) to 
our work, we use the effective source parameters as explained in section A 
of this chapter. The double exponential integral H required for the 
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solution is written in terms of single exponential integrals by Fenlon. 
Merklinger et al. state that Penlon’s expression is an approximation 
accurate to 20^ for Or^SO.l dB, where a is the absorption coefficient 
expressed in dB/unit dlsteuice. The ineqmlity is easily satisfied in our 
experiments. We calculate the single exponential integrals by using a 
polynomial approximation given in Ref. 45. The expansion has a maximum 
error of 2X10”^ for arS4.3 dB. 


C. Theory Applicable to Weak Waves 

A perturbation solution of Burgers* equation for spherical waves 

may be applied to a wave in which no shocks are present. The solution for 

the second harmonic, accurate to second order in €, for a sinusoidal 

signal at r is^^ 
o 


p^r^ a -4a(r-r^) 

o o ^ « o 

^2 “ r 2 ■^22 ® 


(3-10) 


where 


‘22 


■/. 


^ 2[a(r*-r^)] 

S. - - ^ 


-20(r fr 2CO'' 

ar' oe ° / a.-pstoo' 
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It may be shown that, to second order in €, any second harmonic source 
distortion at r^ does not enter the perturbation solution; rather, the 
initial second harmonic distortion diminishes by linear theory as^^ 

/v \ -a (r-r ) 

P 2 (r,t) = pg(r^,t) j e . (3-11) 

The use of instead of allows for the possibility of relaxation 
effects in the medium. Effective source parameters are not required 
because the initial distortion does not affect the nonlinear behavior. In 
general, the phase difference between the two signals represented by 
Eqs. 5-10 and 5-H is nonzero. The total second harmonic is the vector 
sum of the signals. The fundamental component decays according to linear 
theory in this approximation. 


D. Theory Applicable to Strong Waves 

1. Naugol'nyKh et al. (Burgers) Model 

The extra attenuation calculated by Cary using the model of 
Naugol'nykh et al. is applied directly to o\ir work. The extra attenuation 
is 


where 




^/&€k 

sinh ar[l/o^ + ln(r/r^)] 


( 3 - 12 ) 


(3-13) 


An effective source must be defined and applied as in section A of this 
chapter. The solution is valid only in the sawtooth region of propagation. 
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2. Rudnlck Model 

Blackstock has applied the Rudnick model for the fundamental 
to the propagation of a spherical sawtooth wave. For such a wave the 
dimensionless fundamental amplitude (compensated for spherical spreading) 
is 

^1 ~ 1 + a ” 1 + pekr ln(r/r ) 

o ' ' 0 


The nonlinear decay rate is therefore 



(5-l^a) 


while the rate of decay due to ordinary absorption is 



signal 


Adding Eqs. 3“1^ gives the total decay rate 



(3-l4b) 


(5-15) 


Equation 5-15 is a nonlinear ordinary differential equation valid in the 
region r>r. We integrated Eq. 5-15 numerically by the second method of 
Runge-Kutta (see, for example. Ref. 47) . In chapter VI the results sure 
compared with measured data. Effective source parameters are again 
required (see section III.A.). 
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Webster^^ laber found the analytic solution to Eq. 3“l-5 subject 
to the boundary condition of sinusoidal excitation at with a pressure 
amplitude of p^^^. The solution is 


Pi = 7- 


2Pj_o e 


-a(r-r^) 


ar 


1 + pekr e 
o 


[E^(ar^) - Ej^(ar)] 


(3-16) 


where the exponential integral E^(x) is defined (see, for example. 
Ref. 45 ) as 


Ej^(x) = 



The integral must be evaluated numerically. Useful polynomial approxi- 
mptions to the integral are given in Ref. ^5* 



CHAPTER IV 


EXEERIMEOT DESIGN 


A. Introduc'clon 

Certain general considerations in the design of our propagation 
experiment are discussed in this chapter. Of special concern was the 
selection of an acoustic source which could produce measxarable nonlinear 
effects within the physical confines of o\ar vertical propagation path. 

A tower 85 m high, equipped with an elevator which could be used to carry 
a traveling microphone, was available. This tower fixed the length scale 
of the propagation experiment. 

Several sources capable of producing periodic sound of finite 

amplitude in air have been reported in the literature. A few seemed well 

suited for our purposes. To quickly evaluate their effectiveness, we used 

a graphical method for predicting the nonlinear effects associated with 

spherically spreading sound waves. This method has been described else- 

U8,ll 

where but is repeated here for completeness. 

B. Soxirce Selection 

The sources used in our experimental work were chosen on the 
combined bases of (l) operating frequency range and bandwidth, (2) pro- 
jected output levels, (5) materials and equipment already at hand, 

( 4 ) cost, and (5) construction time. The devices considered for possible 
selection are listed below. 

1 . electroacoustic driv»r(s) in three configurations; 
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a. an array of drivers, each with its own small horn 
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5 ^ 


b. a single large horn fed by several drivers" 


: I 


* • 

' t 

-/i 
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c. a dish reflector fed by several drivers ^ 

2. resonance devices in one of the above configurations 

39 ^1 

a. St. Clair generator ^ 

51 

b. stepped-plate transducer 

5. aeroacoustic or air modulated sources 

3T,56 


a. siren 

b. electropneumatic driver 

„ 53 , 5 ^ 
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c. whistle" 

d. Hartmann generator" 

The references cited in the list above point out some o-f the irp^roant 
operating characteristics associated with the scvl '’"s. Our desire for a 
source of easily variable frequency (with perhaps the ability to produce 
narrowband noise for later experiments) quickly eliminated the resonance 
devices, whistle, and Hartmann generator as possible candidates. The need 
to use hi^ audio frequencies, because of their generally favorable effect 
on the development of nonlinear phenomena (see section C), caused us to 
drop the electropneiimatic driver from consideration. The acoustic signal 
from a large horn fed by several electroacoustic drivers may be more 
susceptible to nonlinear distortion within the horn itself than is the 
signal from an eorray of drivers with individual horns. A large horn is 
also difficult to construct, as is a large reflecting dish. We therefore 
concentrated our efforts on the development of a horn driver array. 

From several standpoints, eirrays of horn drivers eu:e attractive 
for use in experiments such as ours. Commercially available drivers, 
which, at least in theory, are reliable uniform devices, may be used. The 


i 
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frequency of operation is easily varied. All electronic control of the 
source is possihle, in contrast to the mechanical apperatus required for 
several of the devices listed above- We were fortunate to have on hand a 
large number of horn drivers cf exceptionally wide bandwidth. Their use 
was attractive not only because of cost but because the wide bandwidth 
feature could be utilized in subsequent experiments on finite amplitude 
noise. 

When it turned oit that the horn arrays generated less intense 
sound than desired, a siren was constructed as an alternate, much stronger 
source- An electrically driven siren as used here may cover a wide range 
of operating frequencies. Sirens with conversion efficiencies approaching 
^ 0 ^ may be built. A siren is less convenient to use than a driver array 
because of the air compressor required. The much higher sound levels 
obtainable from the siren more than compensate, however, for the added 
inconvenience. The horn driver arrays and siren are discussed in detail 
in chapters V and VI, respectively. 

C. Estimate of the Importance of Nonlineeirlty on the Propagation of 

Spherical Waves 

We sou^t to evaluate and maximize the nonlinear effec s 
e. .pected for a given source by specifying design parameters. By 
examining Eqs. 5-5* one recognizes the importance of the two dimensionless 
quantities end These quantities determine the severity of the 

nonlinear behavior. Nonlinear effects may be increased by 

1. raising the source amplitude e, 

2. increasing the frequency f, 

• increasing the effective source radius r^, w>^ich generally 
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depends on frequency as veil as actual source size, 

4 . increasing i.e., switching to a more nonlinear mediTom, euad/or 

5. decreasing the attenuation coefficient a, which depends strongly 
on hianidity and frequency. 

Because several of the factors in this list are coupled, it is not 
isoDediately obvious how to optimize a particular experiment for the 
greatest nonlinear effects. A convenient graphical display of the condi- 
tions required for a certain level of nonlinear behavior may be con- 
structed as follows. 

Consider a sinusoidal wave spreading from a spherical source of 
effective radius r^. For a very weak wave, atmospheric absorption damps 
the wave before shocks can form. Quantitatively, this means that the 
computed quantities r and r^^^ have the following relation to each other; 
r «<r. In a very str<^ng wave, on the other hand, shocks form quickly 

fllfitX 

and the wave travels a significant distance before small signal behavior 
begins, that is, r«r . Thus the condition r=r is a dividing cri- 
i^erion, joarking the approximate threshold of the importance of nonlinear 
effects. If the computed quantities f and r ^ ^ are equal, then more 
serious nonlinear effects are expected to be present. These two condi- 
tions may be written in terms of Pekr^ and ar^ as follows; 


-l/pekr^ 

r=r 2 ar = Bekr e 
meix 0*^0 

-5/pekr 

f=r -> kxr = p 6 kr e 
max 0*^0 




It is convenient here to define a quantity closely related to 3 €kr^. The 



"source frequency level," SFL, is defined as 




SFL . SPLj^ + 20 log^o , (4-2) 

where SPL^^ is the extrapolated source level at 1 m and 
frequency in kilohertz. This level is equivalent to the "scaled source 
level" used by Merklinger et al. For air at 20"C, the relation between 
SFL and 3ekr^ is 

SFL = 20 logj^Q Sekr^ + 16T.2 (dB) 

A plot of Eqs. 4-1 cuts the SFL-ar^ plane into three regions. Nonlinear 
effects for a particular source may then be read as very important, of 
possible importeince, or unimportant for operating points above, between, 
or below the two curves, respectively. The SFL chart in Fig. 4-1 shows 
the operating points for the sources used in our experiments. All oper- 
ating points are measured values. The values of r^ were taken as equal to 
R^/2 for each source, respectively. Judging by the position of point A, 
for example, we expect array No. 1 (seven horn drivers) to be a moderate 
sovjrce in terms of nonlinear effects. Array No. 2 (ten horns) should be a 
stronger source. The siren is indicated to be a very strong soiarce. The 
propagation experiments to be described generally confirm these expecta- 
tions. We return to this subject in chapters V and VI. 


While Rq is commonly applied to piston sotirces, we also use it in this 
thesis with our horn arrays and siren. For consistency, we agree to cal- 
culate Rq on the basis of the "active area" of the source (e.g., the sum 
of the areas of the horn mouths in an array) . 
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FIGURE 4-1 

DESIGN CURVE TO ASSESS THE IMPORTANCE OF 
NONLINEARITY ON SPHERICAL WAVE PROPAGATION 
(CONDITIONS ARE FOR AIR AT 20"C) 

Adopted from Thoebold, Wobitor, and BIreksleek (1976). 
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D. Other Deslm ConslcLeratlons 


Two other Important design considerations were the length scale 
of the experiment euid the effects of source and receiver directivity. Use 
of the SFL cheurt gives a quick indication of the magnitude of the nonlinear 
effects expected hut yields no information about the distances over which 
these effects occur. Once a design for a source appeared promising, 
therefore, it was then uecessary to compute the actual values of r, r, emd 
r^ (see Eqs. 5~5)* Ih was important that some or all of these distances 
be within the limited propagation range available to us. On the other 
hand, an experiment in which all interesting effects occur within a short 
distance (e.g., within 1 m) would not demonstrate whether nonlinear propa- 
gation effects are indeed observable in the innomogeneous medium. 

As for directivity, there were three main considerations. First, 
high source directivity is desirable to boost the SFL (note that otr^ is 
also affected) . Second, the serious effects of refraction by atmospheric 
irregularities on a beam that is too narrow may render the microphone 
measurements very difficult. Finally, although high frequency is 
favorable with respect to SFL, an upper limit is set by the receiving 
microphone, which becomes very directional at frequencies above approxi- 
mately 42 kHz. At higher frequencies, both phase and amplitude response 
are poor. The time waveforms are therefore inaccurate. The sources were 
thus selected and designed to meet the following minimum requirements: 

(l) shock formation should occur within reach of our microphone, (2) tlie 
fundamental frequency of the sources should be in the 8 to 12 kHz range. 



and ( 5 ) the beamwidth should not be less than about 10°. These 
requirements were approximately satisfied by the etrrays and the siren 
that were actually used. 
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CHAPTER 5 
ARRAY EXPERIMEMTS 


A. Introduction 

In this chapter an account of the evolution of the horn arrays 
and the experimental results obtained through their use is given. The 
component testing and selection process is explained. The array transmit 
system and the receive system used for all of the propagation experiments 
are described. Finally, experimental results are compared with theoreti- 
cal predictions. 

B. Array PevelOTment 

A number of commercially available high power horn drivers were 
considered for possible use as array elements. Eventvially three drivers 
were selected for detailed tests: the James B. Lansing, Inc., models 

No. 2470 (rated at 25 W electrical input*) and No. 2482 (60 W) and a 
375“H retrofit ced with an alumintun diaphragm (JO W). Several of these 
were ruled out as being too low in power handling ability, being too 
limited in frequency response, or having excessive internal distortion. 
The frequency response, efficiency, and (unclamped) electrical input im- 
pedance were measured by coupling each driver to a 5 cm i.d., JO m long 
progressive wave tube. The tube was terminated anechoically an i. is 
described in detail in Ref. 6. See Pig. 5-1. The tube's first nonplanar 
mode occurred at approximately 5950 Hz. The receiving microphone was 

-If " ~ 

Manufacturer's ratings are for an rms power input for a continuous 1 kHz 
sine wave. 
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FIGURE 5-1 

PROGRESSIVE WAVE TUBE TEST OF HORN DRIVER 
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15 cm from the driver mouth. The equipment involved in the tests included 
the following. 

1. Hewlett-Packard 5580A spectrum analyzer vith a camera. The 
spectrum smalyzer's sweep oscillator provide ji swept or fixed frequency 
sine wave of constant voltage that tracked with the analyzer center fi-e- 
quency. Measurements were made using a 50 Hz bandwidth. Rated analyzer 
accuracy was ±0.5 dB for the frequency range 20 Hz to 20 kHz. 

2. B&K ^156 1/4 in. pressure microphone and cathode follower. The 
response was rated as flat to within ±1 dB from 45 Hz to 50 kHz, but the 
integration effect caused by flush mounting the micropdione in the tube 
wall moved the hi^ frequency 1 dB point down to approximately 18 kHz.^ 

5. B&K 2108 microphone power supply. 

4. Hewlett-Packard 555D step attenuator. Not shown in Fig. 5“1> the 
attenuator reduced the sweep oscillator signal to avoid overloading the 
power amplifier input. 

5. DHL pulser. This signal gate was activated dviring pulsed mode 
tests of the driver's power handling capabilities. 

6. Altec 25O-B power amplifier. The unit's measured frequency 
response was flat to within ±0.2 dB from 500 Hz to 10 kHz. The maximum 
rms power output was 25O W into a 2.5 0 load. 

7. Pearson model No. 110 current transformer. The current trans- 
former, when fcllowed by a voltmeter, provided ac current data for 
impedance measurements. It also provided eui indication of driver overload. 

8. Hewlett-Packard 400EL ac voltmeter. The voltmeter alone also 
monitored the potential applied for Impedance measurements. 
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All three drivers failed at their rated input powers after 
various periods of operation. Two kinds of damage resulted. Either the 
voice coils physically distorted from overheating and shorted out or the 
fusible link in the voice coil circuit opened. The voice coil was ruined 
in either case. Pulse operation was also tried. The drivers were pulsed 
at approximately double their continuous power rating; some failures 
occurred in this operating mode. The 575"H driver eventually chosen as 
the element for the arrays was derated to 25 Wrms input with fair 
reliability. 

The driver found to have the smoothest frequency response 
(constant voltage drive) in the 5 bo 10 kHz range, the greatest bandwidth, 
and the hif^est efficiency at these frequencies was the 575"H. An added 
advantage of this driver was th.?- availability (from government surplus) of 
approximately 20 units in various stages of disrepair. Fifteen of these 
units were cleaned and refitted with aluminum diaphragm voice coils at the 
James B. Lansing, Inc., factory in Los Angeles. The frequency response 
for a typical 575 “H driver coupled to the -rlane wave tube is shown in 
Fig. 5“2(a) (the response curves in Figs. 5-2(b) and 5”2(c) are discussed 
later in this section) . The test was conducted at low levels. It is 
seen that the response is flat to within dB over the range 6 to 
8.5 kHz €uad is flat to within ±4 dB over the range 0.5 to 8.5 kHz. At 
8 kHz the electric-to-acoustic conversion efficiency (at low levels) was 
found to be approximately 17^. The efficiency varied from unit to unit. 

^ For example, another driver was 13 ^ efficient under the same conditions. 

This difference reflects a 1 dB decrease in SPL. 
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FIGURE 5-2 

FREQUENCY RESPONSE OF JBL 375 HORN DRIVER UNDER VARIOUS CONDITIONS 
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The 5T5“H driver was selected for the propagation experiment. 
An appropriate horn then had no be chosen. A circular exponential horn 
was available from government surplus in sufficient quantities for an 
array. The performance characteristics of the horn required investiga- 
tion. The horn parameters are listed belcnr. 


TABIE V-1 
HORN PARAMETERS 


Flare constant 
Cutoff frequency 
Horn length 


m = 0.0^7 can”^ 

f = 260 Hz 
o 

L = 20 cm 


Horn mouth 
Horn throat 
Driver throat • 


Diameter 


2a = 16 cm 


d = 6.0 
o 

d^ = h.e 

d. 


Area 

A = 190 cmP 
m 

A = 28 
o 

A, = 18 

d. 


The discontinuity in areas between the driver emd horn is ignored in our 
calculations. The farfield acoustic output was later seen to actually 
decrease when a matching section was added to eliminate the discontinuity. 

We based our predictions of the acoustic output of the hom-hom 
driver combination on two related assumptions. First, the sound within 
the horn was taken to propagate as a quasiplane wave. Second, the horn 
mouth was assumed to radiate sound as if it were a uniformly vibrating 
baffled piston. Both assumptions are commonly used in euialyzlng the 
radiation from horns (see, for example. Ref. 55)* Use of the first 
assumption reveals that, when the transmitted wave reaches the horn mouth. 
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its SPL should he 

SPL„ - SPL^ - 20 JIJT^ , 

where SPL^ is the level at the horn throat. If the pressxire due to the 
transmitted wave then causes the layer of air at the horn mouth to vibrate 
as a piston, as implied by the second asstunption, the SPL at a distance r 
in the farfield should be 

SPL„ - SPL„ - 20 106io{^ 

where (=ka^/2) is the Raylei^ distance. At R^ itself, the SPL is 
expected to be 



SPI^ =SPL^-A , 
o 

where A is the difference, expressed in decibels, between the true axial 
pressure at R^ (s. ; Ref. 56 ) and the farfield asymptote. 

Unfortunately, measurements did not confirm the predictions 
based on the two assumptions. Measurements were made at 8 kHz, at which 
frequency AfI dB. If one starts with the value SPL^=15^ dB projected on 
the basis of measxirements in the progressive wave tube, the values 
expected from the formulas above are SPLj^=l44 dB auid SPL^ =1^5 dB. The 
values actually measured were SPLj|j=1^3 ^ =151 dB. Although the 

predicted and measured values of SPL^^ seem close, the good agreement is 
largely fortuitous. The microphone at the horn mouth picked up not only 
the transmitted wave but also the wave diffracted from the mouth rim. The 
strong effect of the diffracted wave was observed: the measured SPL was 
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extremely sensitive to microphone position. In any case, however, the 

discrepancy between the predicted and measured values of SPL^ was 12 dB, 

^0 

a very serious loss. Extrapolating this same loss to an array of seven 
elements operating at 8 kHz, one finds that the value of SFL is reduced 
from an expected value of 172 dB to l60 dB. Point A in Fig. 4-1 corre- 
sponds to the l60 dB value. It can thus be seen that, if the 12 dB loss 
had not occurred, the 7-element array would have qualified as a strong 
source . 

In the hope of diagnosing and curing the cause of the low 
acoustic output, we examined in some detail the axial press\ire distribu- 
tion and the pressure distribution across the mouth of a single horn and 

* 

driver. These experiments were performed in an anechoic chamber. The 
axial pressure distributions for the single horn and for a uniform bafflea 
piston equal in area to the horn mouth are shon-n in Fig. 5-5* Data set 1 
was taken with an l/8 in. microphone (see section D.4.) moimted on an 
optical bench; data set 2 was taken using a l/4 in. microphone mounted on 
the overhead monorail support of the anechoic chamber. The theory and 
data are matched in the farfield because of our uncertainty of the 
behavior of both in the nearfield region. The 12 dB discrepancy between 
predicted and received levels is therefore not evident in the figure. It 
is, however, obvious that the measured nearfield behavior does rot corre- 
spond to the nearfield behavior for a uniform piston. The theoretical 
series of maxima and minima predicted by theory were replaced in practice 
by a monotonic decay in the pressure. Trouble with our use of the uniform 

See Appendix A for the construction and performance details of the 
chamber. 
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FIGURE 5-3 

AXIAL PRESSURE DISTRIBUTION FOR SINGLE HORN AND DRIVER 
(THEORY MATCHED TO DATA IN FARFIELD) 
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piston model is indicated. The Rayleigh distance is, however, confirmed 
as indicating the onset of spherical spreading. 

The liressure distribution across the mouth of the horn is 
compared with the distribution expecteci for a \aniform piston at a 
slightly lower frequency (7.5 kHz) ii. Vig. The asymmetry of the data 

in comparison with theory is glaring. Evidently we do not have uniform 
excitation of the air layer at the horn mouth. The existence of cross- 
modes within the horn is a likely explanation for the observed behavior. 

It is appropriate to mention here three additional aspects of 
radiation from the horn. First, horns of finite length resonate at certain 
frequencies just as cylindrical tubas do. The resonance frequencies may 
be calculated if we may assume an '•.ppropriate output Impedance at the horn 
mouth (see, for example. Refs. 55 and 5 &) . The only assumption considered 
in the literature is the one corresponding to the piston model of the horn 
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mouth. Fo. ■ owing Olson, we write the input acoustic impedeuice at the 
throat of an exponential horn as 


*^ 0*^0 \Zj^tcos(bUe)] + Jp^c^lsin(bL)] 


JA^Z^T^inlbL) 


+ p c lcos(bI 
^o o ' 


where 0=arctan (m/2b), b=l/2 v^k^-m^, and is the acoustic impedance at 

the horn mouth. The driver throat is effectively an extension of the horn 

back to the diaphragm. By examining a defective driver's internal struc- 

2 

t\ire, we found that the throat area at the diaphragm is 11 cm . Ignoring 
the discontinuity b >. the driver horn connection, we plot in Fig. 5“5 the 
input resistance for an equivalent horn, beginning at the diaphragm. The 
resonance peaks predicted at 6470 and 8550 Hz are very close to measured 








FIGURE 5-5 

ACOUSTIC RESISTANCE AT THROAT OF FINITE EXPONENTIAL HORN 


AS-76-n33 
MAT . DR 
n -19-76 



53 




i 



. t 


* , 
i 





5i 

i 1 


ft 


t 

i 

f 


response peaks of the driver arrays. A response peak at 6655 Hz was 
evident when ten drivers were packed together in an array. A peak cen- 
tered at 8250 Hz appe£ired in the frequency response for a single driver 
(see Fig. 5”2(c)) and also for an array of seven drivers. Of course in 
our calculation we have ignored the driver characteristics, which . 
strongly influence the relative magnitudes of the observed peaks. 

The second additional aspect of horn radiation is the possi- 
bility of nonlinear losses within the horn. If quasiplanar propagation 
within the horn is again assimed, it may be shown that the nonlinear loss 
should not exceed 0.1 dB in our case. 

Finally, we must give some comment to the two very different 
frequency response curves exhibited in Fig. 5-2(b) and (c). The presence 
of a difftracted \fave from the mouth rim of the horn in combination with 
the direct wave throu^ the horn allows a pleasing qualitative explanation 
for the characteristics shown. Because the relative phase of the two 
waves within the nearfleld varies widely with the product ka, we sure not 
surprised to see the many peaks and valleys shown in Fig. 5-2(b). The 
relative phase changes sure mch smaller in th "r-rfield sind thus the fre- 
quency response is smoothed (Fig. 5-2 (c). It / be shown (see Appendix B 
or Ref. 56 ) that the direct and diffracted waves combine on the axis of a 
uniform piston to form the time derivative of the source signal. For a 
sinusoidal sovirce, differentiation amo\ints to a multiplication of the 
source signal by Joj, where to is the angular frequency. The result is a 
6 dB/octave Increase in the frequency response. 

It is concluded from the tests and analyses of radiation from 
the horn that, at the high operating frequency we wished to use, 
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crossmodes larobably exist within the horn. Because propagation throuf^ 
the horn is not planar, the excitation of the air layer at the horn mouth 
is not uniform. Therefore the uniform piston model for the horn mouth is 
not reliable. However, the Raylei^ distance still approximately indi- 
cates the onset of spherical spreading. 

'The disappointing performance of the hom-hom driver 
combination touched off a search for an alternative design. Various horn 
modifications and substitutions were tried, including 

(1) shortening the horn to reduce any €uatiresonance effects, 

(2) switching to a conical horn, 

(3) adding a sectoral divider in an attempt to control crossmodes, 

(4) adding a matching section to eliminate the mismatch between 
driver and horn (this modification actually degraded performance by 
2 dB), and 

CQ 

(5) switching to a one-quarter wavelexigth tubular coupler. 

None of these devices yielded more than a 2 dB improvement over the 
original horn. This result was attributed to the higpi frequencies 
involved. At 8 kHz the wavelength, about 4 cm, is small compared to the 
mouth diameter, 16 cm, and to uhe horn length, 20 cm. 

Despite the disappointing performance of the horn- loaded 575“H 
drivers, it was still felt that such an array offered promise as a higjh 
intensity source. The source level could be increased by increasing the 
number of elements. The limit to the size of the eu*ray was set by cost 
and the requirement that the beam not be too narrow. Groups of seven euid 
ten drivers were operated as high intensity arrays. See Fig. 5-6. 

Because of the size of the driver cases, the horns could not be packed 
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with their mouths touching one another. The arrays were therefore 
enrranged to provide the most dense packing of the horns possible. Seven 
drivers (array No. l) were formed in a circular array, 0.5 ni in diameter. 
The beam from this array was expected to be axisymmetric and similar in 
angular dependence to the beam from a circular piston of area equal to 
7 A^. Ten horns were assembled in an equilateral triangle formation, 

0.5 m on a side. The beam from this array (strray No. 2) was expected to 
be asymmetric and somewhat narrower than the beam from array No. 1. The 
gaps between the home in array No. 1 were initially filled by a plane 
circular baffle (l.l m in diameter). The measured farfield SPL on axis 
actually Increased, however, when the baffle was removed. No completely 
satisfactory explanation has been offered for this fact. We speculate 
that a resonance may have been "detuned" by the presence of the baffle. 

In any case, no propagation experiments were conducted with baffled 
arrays. Array No. 1 was operated at 8.2 and 6.6 kSz, while array No. 2 
was operated at 6.6 kHz. These frequencies corresponded to peaks in the 
frequency responses of the arrays. The respective source levels are cited 
in sections F and G of this chapter. 
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C. Array Transmit System 

The electronic equipment comprising the transmit system for the 
arrays is described below. The block diagram appears in Fig. 

1. H^lett-Packard 3^30A spectrum analyzer. The sweep oscillator of 
this device served as an extremely stable variable frequency oscillator. 
This device is labeled "oscillator" in Pig. 5-7. 

2. Hewlett-Pacl»rd 550C step attenuator. 
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FIGURE 5-7 

HORN DRIVER TRANSMIT SYSTEM 
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5. DHL pulser. This signal gate was activated when pulsed 
measxirements were desired. 

4-. Balance controls and input treuisformer. The input transformer 
was needed to prevent the floating ground huss of the amplifiers from 
shorting to earth ground throng other equipment. The balance controls 
allowed individual gain adjustments for each aiiiplifier channel. Thus 
partial compensation for VEuriatlons in the driver sensitivities was avail- 
able. Because the sensitivity variations appeared to decrease when the 

drivers were assembled in an array> this feature went unused. 

60 

5. RCA HC 200H power amplifiers. Each driver was powered by an 
individual Integrated circ\iit amplifier module. Measured total harmonic 
distortion at 1 kHz was 1.4^ for 59 W into 15 

6. Trygon model No. RS^-10 and Harrison 6438 b dc power supplies. 

The supplies provided ±37 *5 V at 5 A to power the amplifiers. 

7. Pearson Electronics model No. 110 current transformer and 
Hewlett-Packard 400 EL ac voltmeter. These two devices were used together 
to monitor the ininit current to the drivers. 

8. Tektronix RM45A oscilloscope. 

D. Receive System 

A block diagram of the receive system used for both the array 
and siren experiments is shown in Fig. 5*8. The receive system x>erformed 
several functions. The absolute harmonic levels and the time waveform 
were measured at each designated propagation distance. The acoustic out- 
put of the source was also monitored against possible equipment drift and 
failure. Beeun patterns of the sound from the various sources were 


recorded. 
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Because of the size of the sources and the anticipated long 
duration of the experimental program, simplicity and stability were 
important factors in the experimenta] arrangement selected for the system. 
As stated earlier, a vertical propagation path was chosen to minimize 
ground reflections. A propagation path parallel to the radio tower 
allowed the use of the tower's elevator as a movable microphone stsmd. A 
single traveling microphone sampled the sound pressure along the axis of 
each sound soiirce. The use of a single measuring microphone avoided any 
response and calibration differences inherent in using several micro- 
phones and simplified the experimental arrangement. The sources were 
ground-mounted for stability. The receive system Included the following 
elements: 

(1) tower with elevator, 

(2) source monitor microphone, 

(5) traveling microphone and anemometer, and 

(4) signal measuring eqtiipment and the enclosure. 

1. Tower 

The tower defined the propagation path by supporting the 
elevator on which the receiving microphone was mounted. See Fig. 5”9* 

The tower was of tubular steel construction and bad s triangular cross 
section with 0.75 ® sides. Guy wires were secured to the tower at mid- 
hei^t euid near the top. The maximum allowable propagation distance 
was 85 m. 

The elevator was an open car of dimensions 0.8 m in width by 
0.9 ® in depth bv 1.2 m in hel^t. The maximum gross load was conserva- 
tively rated at 136 kg. The car was powered by a ground-mounted winch; 



I 



62 


all elevator controls were at gro\ind level. The elevator door was of 
perforated steel. We covered the door with 15 cm thick industrial fiber- 
glass (Owens Corning R-19) to minimize sound reflections to the micro- 
phone. The elevator's position, and therefore the propeigation distance, 
were determined in a crude yet reliable way. A light 100 m long cable was 
calibrated against a measuring tape. This cable was then clamped to the 
elevator and allowed to dangle to the ground where it provided a direct 
measure of the microphone height above the source. Height measurement 
error is estimated as 1^ \inder calm conditions, increasing to perhaps ^1o 
under the strongest winds encountered. The elevator was unfortunately too 
noisy and moved too fast to allow a contiuuous-t racking mode of operation. 

The somd source mount was an indexing table, 91 cm in diameter, 
that was located at the base of the tower. The table could be rotated and 
tilted for adjusting alignment and for measuring beam patterns of the 
sources. The angles of ^be measured beeun patterns were corrected for the 
small displacement of the pivot point below the source. A baffled array 
is shown in place on the Indexing table in Pig. 5-9(13). 

2. Source monitor 

A fixed monitor microphone in. B&K model No. 4136) was 

mounted off-axis in. the farfleld of the source. Hie purpose of the moni- 
tor was to Indicate variations in the acc_stlc output. Depending on the 
source dimensions, the monitor was mounted either on a ring stand that 
stood on the indexing table or on a 1.2 m boom that projected from the 
adjacent utility pole. The utility pole (4.5 m hi^) was covered with 
fiberglass to prevent reflections. The monitor microphone output was 
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checked by switching it into the signal measijring branch equipment just 
before a measurement was made with tne traveling microphone. 

3. Traveling Microphone and Anemometer 

The traveling microphone was mounted on a boom extending 2.6 m 
from the elevator car. The boom was used to remove the microphone from 
the reflection field of the elevator and tower. The microphone-tower 
separation was 3»5 Thus when the microphone was midway up the tower 
the minimun time delay of signals reflected from the tower was 1.7 msec. 

No definitive reflection tests were made. However, if reflections were 
Indeed important in the propagatlo*' experiments, they would have been 
manifested as anomalies in the received cw signals. Because the measured 
time waveforms agree quite well with computer predictions (see Figs. 5 •''6 
and 6-9), the role of reflections in the data was probably negligible. 

The traveling microphone was a l/V in. B&K type No. ^156 or a 
1/8 ir ■- '<■ type No. 4138 mounted in such a way that the sound waves 
impi . ;^:azing (90’) incidence. The microphone was calibrated abso- 

luteo.y to within ±0.2 dB by using a B&K type No. 5220 pistonpbone. The 
output signal was fed to the electronic equipment on the groiund by way of 

61 

J.00 m of B4K microphone caDle. The cable had negligible losses. The 

sections of microplion' cable were clamped together to reduce strain on the 

cable connectors caused by the substantial cable weij^t. Windscreens were 

not used because the ones available to us required leirge freefield cc rec- 
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tlons at high frequencies. 

The microphone oom Included provisions for mounting a snail 
laser to aid in source-receiver alignment. However, the laser available 
to us proved to be too weak for daylight use. Therefore, we j'<ed another 
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alignment scheme which was simple yet effective. Witn tne traveling 
microphone 3 ^ above tae 7-element driver array source, a plumb bob was 
used to bring the source (and indexing table) in line with the microphone. 
The microphone was th'^n sent to a h ght of approximately 76 m. ‘The table 
was tilted to maximize the received sound level. The table was then rota- 
ted 90° and the received sound level was maximized again. This position 
was defined as a 0° tilt with respect to the tower; the table was locked 
ohere. This same positirn -was used for the siren. Acoustic alignment of 
the siren was difficult because its major lobe was so broad. On the other 
hand, because cf the siren* s broad beam, alignment of ohe siren was less 
critical than ±'z was for the array. 

To measure wind speed near the microphone position, we mounted 
an anemometer on the elevator car. This anemometer, an Electric Speed 
Ir'^icator Co. model No. F420C, is a wind-driven dc generator whose output 
pc .ntial is proportional to the wind speed. By means of 92 m of stranded 
cable, which had negligible losses, the output signal was relayed tc the 
ground, where it was connected to a calibrated Ester line -Angus model 
No. AW graphic ammeter. The accuracy of the wind speed measurement was 
ap .::oximately ±0.5 kra/h. The time-averaged (by eye) wind speed was noted 
at each propagation distance. Although more accuracy could have been 
obtained by integrating the time lecord output, it was not deemed worth- 
while to do so. 

The air temperature and relative h\amidlty at groiind level were 
measured with a Bacharw.ch sling psychrometer. Several times when the wind 
was fairly calm, the temperature along the tower was found to be essen- 
tially the same as the tempf'rature at the ground. Isothermal (ambient) 
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conditions were therefore assumed for our data analysis. An electronic 
thermometer (Yellow Springs Instriunents model No. Aj) with remote readout 
vill he added to the elevator in the near future to allow a more ccireful 
measurement of ambient temperatiire . 

We note that knowing temperature and wind speed as a function of 
propagation distance is not exactly equivalent to knowing the wind and 
temperatiare profiles as functions of time. We were forced to assvune that 
the gross profile features were stationary with respect to the time 
required to rxm the propagation experiment. 

4, Signal Measuring Equipment and the Enclosxire 

Figure 5“8 and the equipment list below describe the electronic 
portions of the receive system. 

(l) Two l/h in. B&K 41J6 pressure microphones with cathode followers. 
These comprised the monitor and traveling microphones for the driver array 
experiments. The microphone pressure response was nominally flat from 
50 Hz to 50 WIz (1 dB down points) for grazing (90°) incidence. Because 
the ^156 was a pressxare response microphone but operated under freefield 
conditlors, the following manufacturer's rated freefield corrections were 

61 

added to the microphone response; 

Without protective grid: -0.2 dB (9 to 1^ khzj 

+0.5 dB (19 to 58 kHz) 

With protective grid: +1 to 5 dB (17 to ^ Idiz) 

Spectral measurements in these frequency ranges were corrected accordingly. 
Time waveforms were not corrected, although they could have been when the 
signals were processed by the computer. The justification for not making 
.he corrections is described in section G of this chapter. 
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(2) l/8 in. B&K 4158 pressure microphone. This was the traveling 

microphone used for the later siren experiments. Its shorter rise time 
was useful in obtaining more faithful oscillograms of the sawtooth waves 
encountered. The 4lj8 microphone was rated as flat to within 11 dB from 
20 Hz to 100 kHz. No freefield corrections were required to 70 kHz.^^ 

(5) Two B&K 2108 microphone power supplies. 

(4) Ad Yu 20 dB transistor preamplifier. The preamplifier improved 
the signal-to-noise ratio for the traveling microphone. The measured gain 
was 20±0.1 dB from 5OO Hz to 20C kHz. The phase distortion was negligible 
above 100 Hz. 

(5) B&K 2010 heterodyne analyzer. The frequency domain measurements 
were made with this device. Either overall or bandpass levels were 
selectable. The frequency response (linear) was rated at ±0.2 dB for the 
range 10 Hz to 50 kHz and ±0.5 dB for the range 2 Hz to 200 kHz. Bandwidth 
was selectable from 5-l6 to 1000 Hz in steps cf 10 dB. The variable band- 
width feature was helpful in reducing wind noise and wind-induced fluctua- 
tions. The overall SPL measxirements did, however, include a contribution 
due to wind noise. The sinalyzer also contained much needed averaging 
circuitry. The pressure signal was first convei'ted to its true rms value 
then was RC (resistor-capacitor integration) averaged, and finally was 

log converted for direct meter display in decibels. The level displayed 
was thus the pressure signal average (an arithmetic mean), not the average 
of the pressure levels (a geometric mean). The effective averaging time 
was selectable from 0.1 to 100 sec in 10 dB steps. The analyzer filter 
required manual tuning for each harmonic. The tuning was expedited by 
using the frequency coi.mter to monitor the analyzer's beat frequency 
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oscillator output, which tracks with the filter center frequency. When an 
averaging time of 10 sec could he used, about 5 ™in were required to 
meastire the overall and first three harmonic levels. 

(6) Hewlett-Packard 550OA frequency counter. The frequency co\mter 
fulfilled several needs. It was used to monitor the oscillator frequency 
for the driver array (at which time it was connected to the array trsinsmit 
system) or the siren frequency. It was also used to indicate the analyzer 
bandpass frequency, as described in No. (5). 

(7) Tektronix SM45A oscilloscope with camera. The received signal 
fluctuated because of wind and tirrbulence. The oscilloscope records an 
instantaneous waveform, not a time averaged waveform. Thus there is a 
veiriability among waveforms measured at the same distance. Because the 
wind and turbulence increased with hei^t and because of the cumulative 
effect of propagation through a random medium, the variability also 
increased with height. Therefore the waveforms recorded at the larger 
distances (see Figs. 5-l6 and 6-10) are, in particular, to be regarded only 
as representative. The computer elgorithm requires a measured waveform 
for its input. Because slij^itly different input wa'eforms might lead to 
rather greatly different computed waveforms at long distances, we desired 

a very stable source wa eform. Since little variability was found at 
distances of less than 15 ci, a waveform from this region was acceptable as 
the computer input. 

(8) Low pass (lOO kHz cutoff) filter (constructed by D. A. Webster). 
When the traveling microphone was set at the longer propagation distances 
(above approximately ^5 m), the long microphone cable, acting as a radio 
antenna, picked up signals from a local AM radio transmitter (carrier 


68 


frequency of 1.2 MHz). The tower was within the nearfield of the station 
and also within a stsuiding wave field caused by the Balcones Fault. Radio 
reception was un^^redictable but varied with the time of day and the 
weather conditions. Althou^ the radio signal was very strong at times, 
it did not interfere with the frequency analysis of the acoustic signal. 

It did, however, greatly affect the time waveform displayed on the 
oscilloscope. To retrieve the time waveform for the oscilloscope camera, 
we added the low pass filter at the oscilloscope input whenever it was 
required. The filter effectively blocked the radio frequency carrier but 
faithfully i>assed the acoustic signal. The overall SPL measurements may 
have been affected by the radio frequency component. Because the overall 
SPL measurements were of minor importance and were also affected by wind 
noise, the radio frequency component was not filtered from the input 
signal to the analyzer. 

a. Experimental Accuracy 

We estimate here the accxaracy of the absolute SPL measurements 
to be described. For this analysis, we assume that the air is absolutely 
calm. Several components of equipment were involved, either directly or 
indirectly, in each SPL measurement. If we assume that the measurement 
errors attributed to the various components are independent and that each 
error is normally distributed, the stsundsurd deviation of a measurement is 



where s, , s„, ..., s are the standard deviations associated with the 
X 2 n 

individual components. These standard deviations were the following: 
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, . * 

(s^) pistonphone output: 0.2 d£ 

, V * 

(Sg) polarization potential deviation from 200 Vdc; 0.2 dB 

/ \ * 

(Sj) heterodyne analyzer frequency response: 0.2 dB 

/ \ ♦ 

(sjj^) heterodyne analyzer calibration: 0.2 dB 

(s^) meter reading: 0.05 ^ 

(sg) preamplifier frequency response: 0.1 dB. 

The stemdard deviations associated with the total receive system were 

therefore 0.4 dB. In other words, the probability of weasurement errors 

less than 0.4 dB was 68^t. The probability of errors greater than 0.7 dB 

was 10^ (90^ confidence limits). But we must also concern o\xr selves with 

the total experiment accuracy. That is, still assuming calm air, we must 

also account for deviations in the transmitted acoustic signal. The 

deviation of the source output measured by the monitor microphone was no 

greater than O.05 dB for the driver array and 0.5 dB (5 sec average) for 

the siren experiments. Including these deviations yields a standard 

deviation for the complete system of 0.4 dB for the array experiments and 

0.5 dB for the siren experiments. 

b. Enclosure 

The electronic equipment was housed in a 2.1 x 2.7 x 2.1 m 
(outside dimensions) Craig Systems’ Helicop-Kut. The hut provided storage 
space and equipment protection against Inclement weather. Because of its 
potentially airti^t construction, the hut may be used in the future to 
protect the experimenters from unbearably loud sources. 
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E. Procedure for Propagation Experiments 

This section applies to propagation experiments conducted by 
using either type of sound source. The experimental methods were identi- 
cal. The receive system was first calibrated absolutely using the piston- 
phone. The procedure for making a measurement at a given distance from 
the source was as follows. The elevator was moved to the desired distance. 
The source was activated and the output allowed to stabilize. In the case 
of the siren, stabilization required about 1 min. The driver eurrays 
exhibited no drift following activation. The relative pressure level of 
the source was then measured by the monitor microphone. We then switched 
to the traveling microphone and used the heterodyne analyzer to measure 
the SPLs of the overall, fundamental, second harmonic, and sometimes the 
third hetrmonic signal. For small signal experiments, only the overall and 
fundamental levels were measured. Simultaneously, an oscillogram was 
taken and the wind speed was recorded. This process was repeated as 
quickly as possible at the next distance in order to minimize weather 
(mainly wind speed) variations during the experiment. 

In general, the averaging time required for the spectral 
measurements increased with propagation distance from a minimum of 0.1 sec 
to a maximum of 50 sec. This increase refitcted wind speeds which also 
generally increased with height. But the values of wind speed measxired at 
a fixed propagation distance occasional, j varied by as much as km/h. 

The propagation experiments were conducted under a wide rarige of 
wind, temperature, and humidity conditions. The measured wind speeds 
ranged at the ground from 0 to 4 km/h and, at the greatest propagation 
distances, from to 21 km/h. Several of the experiments took place before 
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the anemometer was obtained. We estimate the vlnd speed range at the 
ground for these experiments as 0 to 24 km/h. The wind speeds at greater 
heists are unknown. The propagation experiments took place within the 
temperature range of 24 to 35^0 and the relative hianidity range 58 to 80^. 

The source was shut down between measvirements to protect the 
experimenters from hi^ level sound. The sources operated at levels 
dangerous to hearing, if the exposure is long, even when both eaxpl\igs and 
eannuff-type hearing protectors are worn. 

A special problem was encountered dtiring the siren experiments. 
The siren's high somd intensity prohibited experimentation during labora- 
tory working hours. During the nijpit, small line voltage fluctixations 
(2 to 5 V) caused by the tower's flashing beacon made the siren output 
unstable. We were thus restricted to taking data between 5:00 and 
6:50 p.m., approximately. For these reasons, no data are included here at 
distances greater than 36 m. Improved wiring will be Installed at the 
tower site to eliminate the line voltage fluctuations at night. 

P. Small Signal Behavior 

Small signal experiments were performed to establish a basis 
against which the finite amplitude propagation experiments could be com- 
pared. 

We first needed to determine whether fluctuations due to wind 
and turbulence would so dominate the measurements as to disqualify our 
homogeneous medium theories. Some of our early propagation experiments 
with the arrays provided an answer to this question. We restrict our 
pr< :it discussion to small signal measurements. It was found that, if a 
S' ciently long averaging time were vised, the data agreed very well with 
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predictions based on spherical spreading and atmospheric attenuation 
alone. Agreement was good even though wind gusts of at least 21 km/h 
occasionally occurred at the greater propagation distances and even though 
the array beams were relatively narrow. Thus, it was not necessary to 
account for turbulent scattering and other random media effects. For this 
reason the terms "small signal" and 'linear theory" have been used here to 
imply the effects of spherical spreading and atmospheric absorption in an 
homof-jeneous medium. Two of the small signal propagation experiments are 
now discussed in some detail. 

One of these experiments is reported in Ref. 11. The source was 
array No. 1. Data were taken over the range 6.1 to 76 m. The operating 
peirameters were f=8.25 kHz, R^=5.2 m (calculated value), SFI#=l60 dB, 
ar^=0.046, and SPL^=l4l.5 dB, where r^=R^/2. These conditions define a 
point below the lower curve on the SPL chart (Fig. ^l). Thus very little, 
if any, extra attenuation would be expected. This expectation is ful- 
filled. The propagation data for the fundamental confirm the linear 
theory curve very well. Some second harmonic distortion was observed, but 
it was not strong enough to draw appreciable energy from the fundamental. 
Similar behavior for the fundamental is shown by another experiment with 
array No. 1, described in section G of this chapter. 

The results of a propagation experiment using array No. 2 are 
shown in Fig. 5-10. The operating parameters were f=6.6 kHz, Rp=5-6 m 
(calculated value), SFL=159 dB, ar ^=0.052, and SPLj^=1^3 dB, where 
r^=R^/2. The maximum wind speed recorded, 21 km/h, occurred at a distance 
of 6l m. Point B on the SFL chart (Fig. 4-l) corresponds to this experi- 
ment. The point is seen to lie in the intermediate region of the SFL 
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chart, yet near the lower limit. One might therefore question whether 
this is in fact a small signal experiment. The solid c\arve in Fig. 5*10 
represents the small signal prediction. It is begun at the 6.1 m measure- 
ment because that is a point definitely ir the farfield. The data fit the 
solid curve quite well out to a distance of about 6l m. There is a devia- 
tion of 1.2 dB at the last measurement point, 75 ra. The values of r, f, 
and shown in the figure were computed by taking r^=R^/2. The amount 

of extra attenuation expected to occur between r^=6.1 m and r=24 m is 
0.9 dB. This loss is not observed. The reason for the discrepancy is 
not known. 

The beam patterns for the arrays and a measurement of the axial 
pressure distribution for array No. 1 are now described. Although these 
experiments were performed at maximum output for the sources, the effects 
of nonlinearity or the data are expected to amount to no more than ±2 dB. 
That is, from the results discussed in section G of this chapter, one 
would predict only minor nonlinear losses for the fundamental. The 
measured beam pattern for array No. 1 at 6.6 kHz is shown in Fig. 5-11* 

The wind speed during the experiment was a steady 6.5 km/h. The beamwidth 
(the angular separation between the -3 dB points) is seen to be 8° at 
6.6 kHz. This figure compares very well with the beamwidth of a uniform 
circular piston with an area equal to the array's active area. The theo- 
retical and measured directivity patterns correspond well only for the 
major lobe. The first side lobe of the array is 10 dB below the major 
lobe. The array was axisymmetric and a symmetric beam pattern is assumed. 

The acoustic pressure on the axis of array No. 1 when it was 
operating at a frequen^'.y of 8.25 kHz is shown in Fig. 5-12. Tlie 
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theoreticr.l pressure distribution for a uniform piston of equal active 
area is also shown. The tneory and data have been ratched in tn-. far- 
field. The Rayleigh distance (5.2 m at this frequency) approximately 
indicated the beginning of the farfield zone. The pressijre distribution 
in the neax'field deviated substantially from the predicted values, muci: 
as it was observed to do for a single horn (see Fig. 5“5)« We attribute 
this deviation to the nonuniformity of excitation at the surface of the 
array. In turn, th«- .iOnuniformity is presumed to be at least partly due 
to the present e of crossmodes in the horns. The uearfield .rregularities 
appear to average 0”t in the farfield, just as they do when the source is 
a single driver (see section B of this chapter). Nonlir..-,ar effects are 
not expected to severely influence the nearfield pressuic distribution 
because c’ the short distances and merely local concentrations of high 
IcVwl sound involved. We did not investigate the question of finite 
amplitude effects in the nearfield more thoroughly because of a lack 
of time. 

The measiired beeim pattern for array No. 2 is shown in Fig. 5-15 • 
The wind speed recorded d\iring this experiment varied from 2.5 to Q.5 km/h. 
The array had a beamwidth of ay.proximately 6° at 6.0 ki^z through the 
planes shown in the figure. The array was not axi symmetric and ohe 
asymmetry is markedly reflected in the beam pattern. The narrow major 
lobe may have made the received signal quite sensitive to variation caused 
by wind convection. 
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G. Finite Amplitude Behavior 

Two array experiments performed at maximum output are described 
in this section. 

1. Array No. 1 

* 

When the source for the propagation experiment was £u*ray No. 1, 

the frequency was 6.6 kHz, the Rayleigh distance (calculated value) was 

2.5 m, and the source level was l44 dE. The operating point on the SFL 

chart, point A in Fig. 4-1 (SFL=l6l dB and ar^=0.0070, where r^=R^/2), 

indicates that moderate nonlinear effects should be expected. Jtoreover, 

because of the values of the range parameters — r=ll m, r =850 m, and 

r =35 m — we expected these effects to be observable well within the 
max 

range of our traveling microphone. These expectations were borne out 
qualitatively. At a distance of 9.1 m, both the measured time waveforms 
and the separation between the SPLs of the fundamental and second harmonic 
indicated the presence of a shock. 

The specific propagation data for the fundamental and the second 
harmonic are shown in Fig. 5-l4. Also shown are predicted curves based 
upon models that are probably too simple for the experiment. The propa- 
gation data for the fundamental show no significant deviation from the 
llneeu* theory curve. (Note that, again, this curve is begun at the 6.1 m 
data point.) In other words, no extra attenuation is evident. Yet, 
because of the evidence of nonlinear effects (for example, the presence of 
shocks at least by 9.1 m and the proximity of the second harmonic curve 
to that of the fundamental), one would expect to observe some extra 

The traveling microphone was a l/4 in. model with protective grid. 
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attenuation. Using ordinary weak shock theory (that is. Ignoring 
atmospheric attenuation), one would predict EXDB values of 0.57 dB at 
6.1 m and O.91 dB at 9*1 The increase of approximately 0.3 d£ is not 
observed in the data. The fact that it is not observed may not be signifi- 
cant, however, because O.J dB is less than the experimental error. 

As for the second harmonic, no analytical prediction that is 
really applicable is available. The second order perturbation solution to 
Burgers* equation may be used for weak waves and ordinary weak shock 
theory may be used for strong waves, but the wave strengths for this 
experiment are in the embeurrassing intermediate region for which analyti- 
cal results are very speurse. Nevertheless, an attempt was made to explain 
the data by applying the perturbation results. It is admitted at the 
outset that these results axe relevant only in the shock free region, 
which for this experiment represents only about the first tenth of the 
propagation path. 

In applying the perturbation results, we started with Eq. 5-10 
and let r^=6.1 mj curve (b), shown in Fig. 5-l^> results. To this curve 
we added the second harmonic signal [curve (a) ] that *^s already present 
because of sotirce distortion at r^j see Eq. 3-H* The second harmonic 
caused by finite amplitude effects is generated in phase with the funda- 
mental. If the source distortion were caused by finite amplitude effects 
within the horns, use of the piston model for the horn mouth shows that 
the second harmonic in the far fie Id would be 90* out of phase with the 
fundamental. Since direct measurement showed that the phase difference 

*Note that "source distortion" here includes distortion that took place in 
the nesnrfield as well as in the horns and the drivers. 
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was approximately 0“, we conclude that either the piston model of the 
radiation is inapplicable or that little distortion occurred in the horns. 
The latter conclusion is consistent with our calculation (section B of 
chapter V) that negligible extra attenuation should be associated with 
propagation in the horns. 

While the soluticxi obtained by combining c\ur\res (a) and (b) 
provides a good fit to the data out to approximately I5 m, it is not 
strictly applicable past a distance of about 9.I m. The second order 
perturbation solution does not accurately represent the second harmonic 
signal after a shock has formed. Bata and the perturbation solution 
diverge beyond r=15 m. The only method of solution available to describe 
the propagation for the reg‘ ■jn past shock formation (for a nonsinusoidal 
source) is the computer algorithm. The algorithm was not used, however, 
to provide predictions for this case. 

A much more suitable test of the perturbation solution was 
provided by the data from the first propagation experiment described in 
section F of this chapter (see Ref. ll). (In that experiment the waves 
were sufficiently weak, so that no shock was formed. ) The perturbation 
results are expected to be applicable over the entire propagation path. 

It turned out that the data for the second harmonic fit the predicted 
curve very well. The predicted curve was based on Eqs. 3-10 and 3-U* 

The experiments for which array No. 1 was the sound source were 
completed before the anemometer was Installed. We estimate that the wind 
speeds at the ground were less than 24 km/h. 
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2. Array No. 2 

When the source for the propagation experiment was eurray No. 2, 
the frequency was 6.6 kHz, the Rayleigh distance (calculated value) was 
5.6 m, and the source level was 1^.5 dB. The operating point on the SFL 
chart, point C in Fig. 4-1 (SFL=l65 dB and CCr^=9-6xlO”^, where r^=R^/2), 
Indicates that moderate nonlinear effects should be expected. The effects 
of nonlinearity should be expected to be more pronounced than those for 
the experiment using array No. 1. 

The pi'opagation data for the fundamental and second and third 
harmonics are shown in Fig. 5-15* Tbe data for the fundamental fall below 
the small signal prediction. The maximum deviation from the small signal 
curve is approximately 1.7 4B. It is seen that the second and third har- 
monics grow relative to the fiandamental before they recede. The minimum 
separation of the fundamental and second harmonic is 5* 4 dB, which value 
occurs at l8 m. The minimum separation of the fundamental and third 
harmonic is 11.4 dB at a distance of 9*1 ®. Because the minimum separa- 
tion of fundamental and second harmonic for a sawtooth wave (the limiting 
form of strong waves) is 6 dB, the measured value of 5*4 dB is a curious 
result. Other data from the vicinity of l8 m confirm this value. Thus, 
random error is discounted as an explanation. The reason for this 
behavior is not clear. Asymmetry of the source wave may explain the 
recorded value. 

Theoretical propagation results for the experiment described 
above were generated by using the computer algorithm. The received wave- 
form at r^=6.1 m served as the source wave. Because this distance is 
outside the Rayleigh distance of the array, the effect of diffraction at 
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the horn mouths was expected to bt excluded fl*om the propagation region 
(see Appendix B). The source waveform was prepared for use in the com- 
puter program, as described in chapter III. The waveform was then fast 
Fourier transformed to allow examination of the relative harmonic levels. 
The levels of the second, third, and fourth haimionics relative to th 
fundamental were - 6 . 9 , -13.8, and -I 8.8 dB, respectively. The freefield 
correction for the microphone was +0.2 dB at the second harmonic and 
-O .5 dB at the third and fourth harmonics; it increased to +1.0 dB at the 
seventh harmonic. Because the energy contained in the harmonics v^bove the 
f\indamental was relatively low and other errors in the experiment were 
comparable or larger than the freefield corrections (especially when wind 
noise is accoimted for), the input waveform to the computer was not 
corrected. However, all other measured harmonic levels were corrected for 
the freefield response. 

The computed propagation curves are shown in Fig. 5”15* The 
general trend in the data is visible in this figure. Our predictions are 
too low for the fundamental and second harmonic and are slijpitly hij^ for 
the third heurmonic. The maximum deviation between experimental and com- 
puted values for each harmonic (at the distance the deviation occurred) 
are as follows; fundamental, +1.8 dB (at 76 m); second harmonic, +1.3 dB 
(at 21 m); and third harmonic, - 0.7 dB (at 76 m). A positive deviation 
indicates that the measiared spectral level is the hi^er of the two. The 
prevalence of underestimated levels would seem to be anomalous. That is, 
overestimations could rather easily be explained by additional losses that 
are not accounted for in the computer algorithm. The deviations were not 
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attributed to the selection of the step size parameters in the algorithm. 
We shall reconsider this discrepancy shortly. 

The minimum computed separation between the fundamental and 
second harmonic was a value of 5*8 This value, which occurred at a 

distance of 12 n, is 0.2 dB less than the value predicted for a sawtooth 
wave. Thus the relative levels of an initially distorted wave may possi- 
bly violate the limiting values applicable to an initially sinusoidal 
wave. An analysis by Webstcr^^ applicable to the preshock region may lend 
support to this suggestion. Webster has shown that, for a source with 
initial second harmonic distortion, the separation of the fundamental and 
second harmonic levels is less than it is in the case of a pirre tone 
source. The amount of reduction depends on the relative phase of the 
fundamental and initial second harmonic. This behavior mi^t carry over 
into the postshock region. In any case, because no definitive accuracy 
test of the computer program has been devised (except for the case of no 
atmospheric absorption), the 0.2 dB difference value described may not be 
significant. 

The computed and measured time waveforms and spectra for 
representative distances are compared in Fig. 5-l6. Note that, because of 
the characteristics of the condenser microphone, the electrical output 
signal is inverted with respect to the acoustic signal. The waveform was 
therefore restoi’ed to its proper apisearance through use of the oscillo- 
scope's Inverting amplifier. Hie plot sizes of the experimental and 
computed waveforms are matched for the source waveform. The remaining 
plot axes are proportional in size to the source wave axes, according to 
the relative gain setting of the oscilloscope. The zero pressure line on 
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the oscillograms does not correspond to a grid mark. This apparent dc 
shift, of the measured waveforms relative to the computed waveforms was 
caused by a maladjustment of the oscilloscope's verticf-l position control. 
For ease of comparison, the experimental and computed results for the 
first three spectral levels are repeated in this figure. Theoretical 
results are ^hown as solid bars; the measurements are indicated by 
crosses. 

It is seen from Fig. ^-l6 that the experimental waveforms out to 
a distance of approximately 6l m are more peaked on the positive half- 
cycle than are the computed waveforms; the negative half-cycles match 
quite closely. The experimental waveform at ?6 m has the same wave shape 
as the computed waveform but has a higher amplitude. The differences in 
amplitude between the experimental and computed time waveforms become less 
dramatic when they are expressed in decibels. These differences increase 
with propagation distance, where the effects of the inhomogeneous medium 
were stronger. Because the oscillograms do not depict time averaged 
waveforms, the oscillograms show more variability tha the spectral 
levels. Finally, because the computed results were all based on one 
initial waveform, small phase or amplitude variations in this waveiorm 
may be carried through the calculations. 

A comparison between 0'’r measured values of the extra 
attenuation and the values predicted using the Merklinger et al. model 
was desired. The amplitude and radius of an equivalent spherical source 
of sinusoidal waves must be knov«n for these computations. The relative 
levels of the first and second harmonics at a distance of 6.1 m indicate 
that a shock Ls already present. (The SPL of the second harmonic is 



90 


;) 

ih 


- 6.9 relative to the SPL of the fundamental. The difference in these 
levels for a newborn shock is -T.9 dB.) But the waveform is also asym- 
metric at this distance. Solutions were computed by assvonlng two extreme 
conditions: (l) a sine wave exists at r^=6.1 m and (2) a newborn shock 

wave exists at rj^=6.1 m. (This assumption mean*! that a sine wave exists 
at r^=l.^ m.) These two assumptions lead to the two curves shown in 
Fig. 5-17. The measured data are effectively bracketed by the two predic- 
tions. This result is encouraging, especially when it is recalled that 
the experimental error for calm air is of substantial size on the scale of 
the figure. 

It is appropriate, though somewhat difficult, to determine 
whether the beam from array No. 2 was affected by diffraction. (See the 
discussion in section B.J. of chapter II.) In current theoretical treat- 
ment the amplitude variation across the beam is assumed to be 

Gaussian. Althou^ the variation is not Gaussian for the beam from the 
array, we may very rou^ly approximate it to be so by defining the effec- 
tive beam radius b (see Eq. 2-12) in terms of the points on the measured 
beam pattern (see Pig. 5"15) at which the respr-.je is down 8.7 dB. 

Because the measured beam is narrower on one side than the other, the 
estimate for N, at 1^ m, is a range of values rather than a specific 
value, namely 


0.1<N<0.»V 

Because the beam spreads spherically, b Increases and e decreases with 
propagation distance. The net effect on N (see Eq, 2-lj) is to decrease 
as l/r. 


i. 




II 


i 


\ 


f: 

i 








92 


Several characteristics of the experimental waveforms for our 
propagation experiments axe in qxialitative agreement with the predictions 
regarding the effects of diffraction. The asymmetry of the waveforms 
is of the type predicted to occur, that is, the peaks are shenrpened and 
the trouj^s are rounded. The asymmetric waveform sketched in Fig. 2-1 
corresponds to the case of N=0.4, the upper limit on N estimated for oijr 
experiment. (Note that Pig. 2-1 Is a spatial waveform while the oscillo- 
grams enre time waveforms.) Although the observed asymmetry might be 
attributed to the fact that the initial waveform (at r^=6.1 m) is asym- 
metric, the computed waveforms, which are not affected by diffraction, 
show the asymmetry dying out monotonically with distance. This reversion 
to a symmetric waveform is probably due to ordinary nonlinear distortion. 
The fact that the measured waveforms Increase in asymmetry with distance, 
at least out to a distance of 12 m, implies that an effect not included in 
the computer algorithm is at work. The fact that the asymmetry eventually 
recedes with distance in the measured waveforms is consistent with the 
decrease in the parameter N with distance. When N decreases, ordinary 
absorption and nonlinear distortion become domlneuit, and the asymmetry is 
gradually "damped." When the siren, for which the value of N was much 
less, was used later, the asymmetry did not follow the same pattera. 

(See Fig. 6-10.) Although some initial asymmetry was present, it appeared 
to decay monotonically with distance in both the measured and computed 
waveforms. For these reasons, it is felt that the effects of diffraction 
may have been observed in the hij^ amplitude propagation experiment in 
which array No. 2 was used. 
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We briefly suimiiarize the general i'esults of the array 
experiments. The expectations of wave behavior based on the position of 
the operating point on the SFL chart were borne out. Where theoretical 
predictions were applicable to the experiments, agreement with the data 
was reasonably good. This evaluation applies to the propagation curves, 
time waveforms, and the extra attenuation. The predicted hetnnonic levels 
were generally hi^er than the experimental levels; the measured extra 
attenuation was less than the expected values. Because the predictions 
used are Intended for homogeneous media, it would appear that nonlinear 
effects are not completely neutralized by the random medium, at least for 
the conditions encoxmtered during our experiments. Finally, the asym- 
metric distortion predicted by Rudenko et al. is demonstrated qualita- 
tively in otir recorded waveforms. The details of the effects of 
diffraction on initially asymmetric waves are unknown at the present time. 



CHAPTER VI 


SIREN EKPERIMEIKFS 


A. Siren 

A much more poverful acoustic source was needed in order to 
demonstrate much stronger nonlinear effects than those obtained with the 
arrays. Electrically driven sirens with external air compressors have 
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been used to generate extremely hl^ sound levels in open air.^ ' 

A siren basically consists of two perforated plates, a fixed stator and a 

spinning rotor, suid a supply of compressed air. When the holes (or ports) 

in the stator and rotor align, a burst of compressed air is emitted. 

Because the i»rt disuneter is nozmially much smaller than a wavelength, the 

port behaves like an acoustic monopole. The ports are commonly spaced 

equally in a cix’cle on the rotor and stator plates. The siren is then 
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expected to act as a ring. The frequency with which the ports come into 

alignment, and thus the acoustic frequency, are determined by the motor 

speed f_id the number of plate perforations. The sound pressure output is 

a complex function of the port shape, size, and the pressure difference 
37 

across the plates. 

We elected to construct a siren scaled down in t ze from the 

38 

model described by Allen and Rudnlck.^ The main design criteria were 
that ^:le siren operate at a frequency of approximately 12 kHz emd produce 
i^earby an SPL of I50 to I60 dB. The final output level was to be deter- 
mined by the air compressor available. The mechanical design was by 

The threshold of pain for humans is an SPL of approximately 120 dB. An gg 
SPL cf 160 dB is fo\ind in the nearfleld of a 10,000 lb thrust Jet engine. 
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J. Byers and L. Guyer of Applied Research Laboratories. A schematic 
dleigram of the siren (sketch by W. N. Cobb) is shewn in Fig. 6-1. The 
drive motor was a l/2 hp tool post grinder (Precise model No. Super 6o), 
which had ein vinloaded spindle speed of 45,000 rpm. By mounting the motor 
directly Inside the air chamber, we hoped to use the siren air flow for 
extra cooling of the motor. It was believed that this eirrangement would 
allow the motor to operate above its normally rated load. The rotor was 
machined from tiueuiium for strength and safety. Twenty conical ports were 
drilled in the rotor and stator. The ports were equally spaced with 
centers on a circle of radius q=7*6 cm. The port diameters were: rotor, 

0.87 cm input and 0.53 cm output; and stator, 0.55 cm input and 1.1 cm 
output. Assuming a loaded motor speed equal to 8o^ of its unloaded speed, 
we expected the maximum frequency to be 12 kHz. The siren output was 
coupled to the atmosphere through a 6.4 cm long conical horn (27" flare 
emgle). The completed siren is shown in Fig. 6-2. 

B. Siren Transmi t System 

The siren transmit system, shown in Fig. 6-3, had the followliig 
major compionents: 

(1) General Radio 200 CU Variac. The siren speed was controlled by 
veurying the ac voltage applied to the motor. The upper limit to the speed 
was set by the maximum motor current, which was 5 <6 A. 

(2) Pearson Electronics model No. 110 current transformer and 
Hewlett-Packard 400 EL ac voltmeter. These devices monitored the motor 
c\irrent so that overloading (and overheating) could be avoided. 

(5) Davey 210 WDS-TS piston type air compressor (truck-mounted). 

The air supply for the siren was rated at 94 liters/sec at 6.8 atm (gauge). 
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The supply included an air cooler and storage tank. The compressor could 
supply only enough air to charge the siren chamber to a m avi mu m pressure 
of 1 atm (gauge). Because of pressure regulation problems, the chamber 
pressure could be reduced only to approximately 0.5 atm. 

C. Siren Operation 

Several problems were encountered during the initial testing of 
the siren. The loading on the motor was much higher than had been expec- 
ted. We assume that the loading, which was found to increase with chamber 
pressure, was caused by air drag on the rotor. The design frequency, 

12 kHz, could not be achieved. The meiximum (stable) frequency was about 
6 kHz at full chamber pressure and 8 kHz at minimum chamber pressure. At 
maximum acoustic output (full chamber pressmre), the SPL of the funda- 
mental (6.5 kHz) was 165 dB at a distance I9 cm from the stator plate. 

More seiious than the drop in operating frequency, however, was the over- 
heating of the motor that occurred when full chamber pressvire was used. 
Operation at a motor current of 6 to 6.5 A (the rated maximum current was 
5.6 a) eventually caused the motor windings to fall. After the motor was 
repaired, the motor current was carefully monitored. With the chamber 
pressure reduced to 0.7 atm, the current was eui acceptable 4^.9 A. Under 
these conditions the SPL (again at 6.5 kHz) at 19 cm was only O.5 dB less 
than the value achieved when full chamber preasijre was used. Thereafter 
the siren operated very reliably. 

At low chamber pressures, which were needed for the small signal 
experlioents, regulation of the air supply became a problem. Variations in 
the chamber pressure caused corresponding variations in the acoustic 
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output. The output variations were reduced to 0.2 to 0.3 dB by diunping 
excess air pressure into the atmosphere euid running the compressor con- 
tinuously. In this way stable operation at a chamber pre.'>sure as low as 
0.3 atm was achieved. 


D. Evaluation of the Siren as a Sound Source 

In this section, we describe the Important acoustical 
cheu:acteristics of the siren. The Rayleigh distance (based on the area) 
of the horn mouth was calculated to be 38 cm for a frequency of 6.1 kHz. 
Although the nearfield was not probed in detail, eno^igh measurements were 
made to show that the waves were already spherically spreading at a 
distance r=19 cm. When the siren vas operated at high amplitudes, a wave 
at this distance was approximately triangular in shape. The wave dis- 
torted into an asymmetric sawtooth by the distance r=38 cm. Finite 
amplitude effects are believed responsible for this distortion. 

Microphone measurements made near the horn mouth were disturbed 

by air flow. We did not attempt to determine whether the dc air flow 

through the siren or the "sonic wind” effect described by Allen and 

Rudnick^® was the main contributor to the disturbance. It is suggested 

that the atmosi^eric wind, which blew crosswise to the acoustic beam, may 
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have acted to disperse emy sonic wind present. (Allen''^^ employed an 
electric fan for this purpose for his indoor experiments.) In any case, 
measurements made at r =38 cm and beyond were seemingly undisturbed by the 
air flow. 


The conical horn attached to the siren provided a gain of 
approximately 8 dB over the siren alone. The gain measurement was made at 
a distance of 19 cm. 
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The beam pattern for the fundamental measured at a distance of 
15 m Is shown in Fig. 6-4. The siren was operated at a high amplitude 
(an of approximately l49 dB) for these measurements. A beamwidth 

of approximately 50° is Indicated. The siren alone should be expected to 
act as a ring source. Data described by Allen and Rudnick for their siren 
are in good agreement with the predicted beam pattern for a ring source. 
Their measurements were made at a short distance (25 cm) . The measured 
beamwidth of their siren increased when an exponential horn was added. 

With the addition of a horn, the source might be expected to behave as a 
baffled piston. The predicted beeun patterns for both uniform ring and 
piston models for our siren are given in Fig. 6-4. The measured besun 
pattern does not confirm either model. However, in ovir discussion we have 
thus far neglected nonlinear effects. At a distance of 15 m, where the 
measurements were made, a significant loss due to nonlinear effects had 
already occurred (see section E of this chapter). It is suggested that 
the major lobe pattern has the blunted appearance characteristic of high 
eunplitude sources. Indeed a 6 to T dB reduction, due to extra attenuation, 
at the center of the major lobe would cause the beam pattern for a piston- 
like source to resemble oxxr measured results. 

An estimate of the acoustic output of the siren at the 
fundamental frequency foi high amplitude oxjeration was desired. The power 
output may be estimated by Integrating the received sound pressvire of the 
fundamental over the wea of the beam. Carrying out this caJ.culation by 
using the information given in Fig. 6-4, one obtains a power of 60 W at 
the fundamental frequency. But, as noted above, extra attenuation has 
probably reduced the level by 6 to 7 dB, at least in the center of the 
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major lobe. If one instead extrapolates the pressure measured on-axis at 
0.58 m (the Rayleigh distance of the horn) oul co 15 and uses the same 
heeun pattern measrared there, a power rating of 4l0 W is o' lined. The 
corresponding extrapolated SPL on-axis is 125.7 dB versus the jl7.5 dB 
value measured at I5 m. Because nonlineeir losses tend to square off the 
major lobe, the value 4l0 W may be an overestimate. If our estimate is 
correct, however, 85^^ (350 W) of the siren’s power output (at the funda- 
mental frequency) was lost within the first 15 m of propagation. In the 
future, for purposes of comparison, a small signal beam pattern for the 
fundamental should also be 'pleasured. 

We were w*. '.hout an airflow me'.<;r during these experiments and 
the siren conversion efficiency was not meat ured. Efficiency was not a 
major concern in this application. 

Because the acoustic output of the siren was an asymmetrical 
sawtooth-type wave even at very short propagation disteinces, we desired a 
check on whether low level sound from the siren would exhibit small signal 
behavior. When the chamber pressure was reduced to approximately 0.3 atm, 
we were able to conduct an experiment at 6.1 kllz at a source level 
SPLjj^:-l*lO dl'. The operating point on the SPL graph (Fig. 4-i) is point D 
(SFL=156 dB, ar^=8.5AlO”^, where r^=R^/2). The pi'opagation curve (see 
Fig. 6-5) for the fundamental matches the linear theory extremely well. 

The deviation is less than 0.9 dB out to a distance of 46 m. 

E. Propagation Experiments at Finite Amplitude 

'Die data from three propo,t',..tion experiments with the siren 


operating at high amplitudes are shown in Figs. 6 6, 6-7 » and 6-8. The 
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valxie of SPLj^ for all three runs was approximately 1U9 cB, where 
fundamental firequencies of 6.1 emd 6.5 kHz were used. The corresponding 
points on the SFL chart are shown in Fig. 4-1. The variation in position 
of the three points results from the vide ranging values of a associated 
with varying weather conditions as well as from the silently different 
fundamental frequencies used. The measiu'ed operating parameters and the 
calculated range x>arameters are shown in the figures. All range param- 
eters are calculated using the effective source parameters, as described 
in chapter HI. The linear theory propagation cinrve for the fundamental 
is drawn from the first data point at or beyond for each experiment. 
Measurements at propeigatlon distances were not considered in the 
data reduction. 

Several important differences between the results of the three 
hi^ amplitude experiments must be noted. First, the siren output during 
experiment 1 was unstable beca\ise of line voltage fluctuations. Therefore, 
no theoretical predictions are calcvilated for this run. Ser'ond, the 
traveling microphone in experiments 1 and 2 was a B&K type 4-156 (l/4 in.) 
with the protective grid in place. A B&K type 4138 (1/8 in.) microphone 
without a grid was used in experiment 5« The euigle of incidence was 90° 
in all cases. Waveforms received at similar distances for experiments 
1, 2, €uid 5 are shown in Fig. 6-9. Althou^ the dlstemces are not all the 
same, one would expect the propagating wave to have about the same wave- 
form in each case. Yet the three waveforms appeeur quite different. 
Analyzing the harmonic components of the waveform from experiment 1 on the 
Hewlett-Packard digitizer, we found that the fifth and sixth harmonics 
(30.3 and 36.6 kHz, respectively) were approximately 3 ^ above the 
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relative levels expected for a perfect sawtooth wave. The frequencies of 
these harmonics lie in a region in which there is a 5 dB response peak of 
the microphone (for 90“ incidence with the protective grid in place). A 
few selected measurements were performed to show that the "enhancement" 
of the harmonics decreased when the protective grid was removed. But to 
ensure that the rise time of the shocks was not severely limited hy the 
frequency response of the traveling microphone, we used the l/8 in. model 
(without the protective grid) in experiment J. It is clear that use of 
the l/8 in. microphone gave much improved results. The "biomps" appearing 
in the output of the l/k in. microidione may probably be attributed to 
diffraction and reflection effects involving the slotted protective grid. 

Third, a most curious result is provided by experiment 2. The 
second harmonic appears to have approached within 3 dB of the fundamental. 
For a perfect sawtooth wave, the second harmonic is 6 dB below the funda- 
mental. We were perplexed by this discrepancy at the time the experiment 
was performed emd sought to pinpoint the error source. The traveling 
microphone, a l/4 in. model with grid, was replaced by a duplicate, but 
the same results were obtained. The meter readings were triple checked. 
Later, analysis of the oscillograms by the digitizer and calculator indi- 
cated that the difference in harmonic levels was actually close to 6 dB. 
We concluded that the vinexpectedly close proximity of the second harmonic 
to the fundamental was caused by a malfunction in the heterodyne an«'' ^er. 
Note, however, that the measured fundamental levels appear to be unaffec- 
ted by the assumed malfunction, xo is also noteworthy that the funda- 
mental second harmonic separation in the other two experiments is close 
to 6 dB. 
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To obtain theoretical predictions to compeire with data from 
experiment 3# we used the computer algorithm discussed in section B.l. of 
chapter III. The waveform measinred at a distance of O.38 m, where signal 
fluctuations caused by the dc air flow were small> was used as the input 
for the computer program. The freefleld corrections for the traveling 
microphone are less than O.5 dB out to a frequency of 70 kHz. The input 
wave was therefore not ^reefield corrected. The sizes for the distortion 
steps and the absorption correction steps in the computer program are 
constant. Because the propagation distance is nearly two orders of mag- 
nitude greater than the initial distance, a considerable amount of 
computer time was required for program execution. In computing the 
results for the first five range points, we used an absorption correction 
step size of 8 cm. The program was executed a second time with a step 
size of 61 cm to calculate predictions for the last three points. This 
change of step size apparently caused an offset in the computed results 
(see Fig. 6-8). The problem may be remedied in the future by allowing the 
step size to increase exponentially with propagation dlstemce. 

The computed curves and experimental propagation data for 
experiment 3 are compared in Pig. 6-8. (The data and predictions are 
also displayed in a different fashion in the third coliamn of Fig. 6-10.) 
Agreement between theory and experiment appears to be qxiite good, although 
there is some indication that the agreement may deteriorate at large 
distances. The change in step size described in the previous paragra|di is 
indicated by the break lines in the figure. The maximum deviations 
between computed and experimental results are as follows: fundamental, 

-1.6 dB (at 30 m); second harmonic, -1.4 dB (at 30 m); emd third harmonic. 
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-3*7 ^ (at ^0 m). The negative deviations indicate that the measured 
values are lower than the cosrputer predictions. The measurements at JO m 
may have heen affected by line voltage fluctuations caused by the auto- 
matic activation of the tower beacon. 

The minimum separation of the fuztdamental and secoad harmonic 
measured in experiment J vas J.J dB (at 6.1 m). This value is again 
within approximately one standard deviation of the value expected for a 
sawtooth wave. The computer prediction at 6.1 m shews a separation of 
J.4 dB. Thus it may be possible for the fundamental- second harmonic 
separation to drop below to 6 d£ for an initially aisyniBetrlc stixaig wave. 

The separation between the linear theory curve and the data for 
the fundamental represents the extra attenuation. Values of extra attenu- 
ation of up to 7 dB are evident in Fig. 6-8. An interesting comparison 
between experiment J and the small signal run may be made. Althou(^ the 
source levels for the two experiments are about 6 dB apart, the received 
levels at JO m are within 1 dB of each other. The question arises as to 
Aether the condition of acoustic saturation is approached in this experi- 
ment. Before considering the subject of saturation, however, we first 
describe the computer predictions for the time waveforms in experiment J. 

The experimental and computed time waveforms and spectra for 
representative distances for experiment J cuo> compared in Fig. 6-10. 
Agreement between theory and experiment appears to be good, particularly 
at the shorter propagation distances. Again recall that Instantaneous, 
not average, time waveforms are compared. 

The value of N, the parameter gauging the relative Importance of 
beam diffraction to nonlinear distortion, is estimated to be very small 
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for this experiment^ namely NsO.023 nt 1^ m. I^erefore, little, if any, 
beam diffraction effects are predicted. Although asynnetry is present in 
the waveforms, it is attributed to Inherent characteristics of the siren 
source. Note, for example, the fact that the asymmetry is definitely 
greatest close to the source. Furthermore, in contrast to the asymaetry 
of the waveforms in the array No. 2 experiment, the asymmetry here decays 
rapidly and monotonlcally with distance. 

We now take up the matter of acoustic saturation as related to 
our siren experiment. The maximum amplitude attainable at a given 
distance, regardless of the acoustic power available at the source, is 
termed the saturation amplitude. While saturation theory for plane waves 
has met with considerable success, the theory for si^erical waves has 
lagged behind. However, Laird (see, for example, Ref. 2) has found the 
saturation amplitude for diverging sj^rical waves in the sawtooth region. 
Bis result, adapted for the fundamental component of the acoustic pressure, 
is 

a 

As Shooter et al.^ have pointed out, this formula is not expected to be 
accurate in the region beyond Shooter et al. proposed an improved 

formulation which takes small signal absorption into account for the range 
beyond r . However, because weak siiock theory affords a reasonably 
accurate description of the received level of the fundamental for our 
experiments, even at distances beyond (see below), we calculate the 
saturation level of the fundamental as a function of r on the basis of 
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Eq. 6-1. Tlie value of is determined by the extrapolation procedure 
outlined in section A of chapter III. The results are plotted iti 
Figs. 6-11 and 6-12 for experiments 2 and 5> respectively. The measured 
data approach within 1 dB of the predicted saturation levels. By impli- 
cation, therefore, an increase in acoustic output near the siren would 
hardly be noticed at distances exceeding, say, 5 m. The bounding of o'jr 
experimental data by the saturation level is Indeed an important result. 

The measured level of the fundamental is compared to predictions 
based on linear theory, weak shock theory, and the Rudnick model in 
Figs. 6-11 and 6-12 for experiments 2 and 5> respectively. Both weak 
shock theory and the Rudnick model provide substantially better predic- 
tions than linear theory. The small effect of atmos]^ieric absorption in 
these experiments may be seen by comparing the Rudnick model and weak 
shock theory curves. The difference in the predicted levels at a distance 
of 50 ™ is about 2 dB. Although on theoretical grounds the Rudnick 
model would seem to be preferable, the data do not clearly indicate that 
it is. Additional experiments a<i longer propagation distances and/cr 
higher frequencies are clesurly needed. The effective source pai'ameters 
used in making the calculations for Figs. 6-11 and 6-12 were extrapolated 
from waveforms at 0.55 and O.69 m, respectively. These distances were the 
smallest at which sawtooth waves were indicated. The resulting values of 


^eff* acoustic Mach number of a spherical sinusoidal source 0.' radius 

Ty ^ , are very high. The values are in fact above the amplitude limit 

7 

for the applicability of weak shock theory to plane waves. But, because 


of spherical spreading and the small values of r. , the amplitude Is 

°eff 

quickly reduced to values within the region of applicability. It must be 
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recognized that our extrapolation procedvire was applied froo a farfleld 
point, where a convenient valtie of a could be determined. 'Riere Is some 
arbitrariness Involved In the selection of the value of r^. Because waves 
of very high amplitude quickly become sawtooth waves, however, much of the 
asymmetry In wave shape present neeur the siren Is lost at larger distances. 
Therefore, basing the extrapolation on a farther rather than a nearer 
range point should be expected to Improve the reliability of the 
parameters. 

The measured extra attenuation Is compeured to that calculated by 
the Merkllnger et al. , Cary, and Rudnlck models In Figs. 6-1? and 6-l4 for 
experiments 2 and ?, respectively. The curves for the Cary model end at 
the limiting distance for the validity of the model. This limit Is a 
dlstemce analogous to r ^^ , occurring where 

The extra attenuation Is assumed to reach a plateau value at the distance. 
The reader Is reminded that the extra attenuation predictions plotted in 
Pigs. 6-13 and 6-l4 are referred to the level of the fundamental at the 
distance r^^ from which the effective source parameters were extrapolated. 
The extra attenmtlon curve for the Rudnlck model in Fig. 6-I3 Is an 
estimate only. The Rudnlck model is applicable only In the sawtooth 
region (r20.?3 m In experiment 2). The source penrameters for use with the 
Rudnlck model were therefore extrapolates from rj^=0.55 Parameters for 
the other models were extrapolated from r =0.38 m (Just beyond the R^). 

In order to maintain consistent sets of soiree parameters for all three 
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models and yet reference the results to the level at the range r=0.38 m, 
the JfiXDB predicted hy use of the Rudnlck model vas added to the EXIB 
measured at r=0.^ m. The agreement between measurements and theories Is 
fair to good. There Is a large amount of scatter In the data. The 
experimental error is quite large on the scales of the figures, Just as 
It was in Fig. 3~17> However, all three models are seen to fit the 
observed extra attenuation data to within 2 dB in Fig. 6-13 ^d 1 dB in 
Fig. 6-l4, 





CHAPTER VII 


CONCLUSIONS 

A summary of the results of our e^qperiments is given in this 
chapter. Several further studies are planned for the future, and the 
direction of these studies is indicated. 

We have carried out a series of experiments on the propagation 
of finite amplitude sound outdoors. The main purpose of the eiqperiments 
was to determine the extent to which the outdoor environment, mainly the 
landom Inhomogenelty of the medium, affects finite amplitude propagation. 
Second6ury purposes were to record eind, if possible, e:^laln any new non- 
linear phenomena. The soxirce frequency was generally in the range 6 to 
8 IsHz, and the source level SPL^ varied from l^K) dB to 1^ dB. Two 
different sources were used: an array of either seven or ten horns and 

a siren. In terms of nonlinear effects, weak to moderate waves were 
generated by the arrays; strong waves were generated by ^he siren. The 
propagation path was vertical amd pEirallel to an 8^ m tower, whose eleva- 
tor carried the traveling microphone. Maximum propagation distances used 
varied from 30 m to 76 m. Meterologlcal conditions were as follows: 
temperature range 2k to 35*0, relative humidity range 38 kQft>, and wind 
speed (at ground level) 0 to 24 km/h. The mesisurements v^re made both 
during the day and at night and were mednly done during the months of 
Avigust through October 1978 • 

Several theoretical approaches were used to explain the data. 

A graphical method, the SFL chart, was used to gauge the magnitude of the 
nonlinear behavior esipected. The SHi chart proved to be a quite reliable 
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guide in predicting the gross importance of nonlinearity. Pestorius's 
computer model^ for the propagation of finite aaplitude plane waves was 
modified to include the effects of spherical spreading and atmospheric 
ahsozption. This algorithm enabled predictions of both time waveforms 
and spectra to be made over the propagation path. The ad hoc propagation 
models of Merklinger et al.,^*^ Rudnick-Webster,^^ and Cary^^ were used 
where applicable. Predictions baused on ordinary weak shock theory werc~” 
checked against the siren data. A perturbation solution of Burgers' 
equation was con^ared with data for a horn array. 

It was found that nonlinear propagation effects were modified 
by inhomogeneity of the mediiim in the following way. The atmospheric 
inhomogeneities did cause significant fluctuations in the instantaneous 
acoustic signal. These fluctuations were most apparent at large propaga- 
tion distsinces. The waveforms varied both in amplitude and wave shape. 
However, with sufficient time averaging when the spectral measurements 
were made, the frequency domain data (i.e., the levels of the first three 
harmonic components) largely conflxmed predictions based on theory for 
homogeneous media. There was a slight indication that these predictions 
may be less reliable at the longer propagation distances. Thus additional 
e:q>erlments need to be performed at larger distances. 

The siren-generated sound waves were very Intense and were 
greatly affected by nonlinearity. Intense sawtooth waves were produced 
very near the source. The propagation data for the fundamental came 
within 1 dB of the predicted saturation level (the prediction was based 
on weak shock throry). Extra attenuation up to about 8 dB was measured. 
Values of the extra attenuation (EXDB) predicted by use of the Budnick, 



125 


Merklinger et al., and Cary models agreed with the measured values to 
within 2 dB. The Rudnick model and weak shock theory predictions of the 
level of the fundamental were within 1 €uad 2 dB, respectively, of the 
data for maximum propagation distances of ^6 m. The conqputer algorithm 
predictions of the fundamental component were within 2 d£ over the same 
distances . 

The beam pattern for the siren, meastired at finite amplitudes 
at a distance of 15 m, appeared to have the hlimted major lobe charac- 
teristic of high amplitude sources. The siren's acoustic power output at 
the fundamental frequency was estimated to be 4l0 W. A small signal beam 
pattern for the siren should be measured as a check on these leist two 
conclusions . 

The results of the array experiments were less clear cut. 

The measured beamwidths of the arrays were from 6 to 8*. Nonlinear propaga- 
tion distortion was produced. The maximum value of EXDB measured was 
approximately 1.5 dB, or approximately 1 dB less than the value predicted 
using the computer and Merklinger et al. models. 

The received waveforms for the e:q)eriment in which array No. 2 
was the source were initially asymmetric. Tbe asymmetry first grew with 
propagation distance and then decreased at larger distances. This behavior 
is qualitatively compatible with the speculation that narrowbeam-type 
diffraction was important in the e:^eriment. The computer algorithm, 
which did not contain the effects of diffraction, predicted a monotonic 
decrease of the asymmetry. 

The results of an e^qperiment described in Appendix B show the 
interaction of diffraction and nonlinear propagatior distortion for the 
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radiation of intense sound from a circular horn. Both e3q>eriii>ental 
results and theoretical predictions indicate that, for the conditions of 
the experiment, the diffraction cmd distortion stages of the propeigation 
may be treated separately. Additional studies on this subject are needed 
to e:^lore the cases where the two effects may not be handled separately. 

Our project was planned as the prototype e^riment for 
future studies on finite atiqplitude noise. The extension of this work 
is clear. Many eispects of the propagation of finite amplitude noise are 
experimentally and theoretically uneiqplored. One subject of significant 
practical interest is the change in spectral content of noise with 
propagation distance. This change gives rise to extra attenuation in 
some spectral bands but to apparent amplification in others. Another 
subject of interest is the effect of ground reflection on propEigation of 
finite -aiig>lltude noise. Finally, it would be of Interest to carry out 
a basic study of the local effects of nonlinearity in the Intense sound 
field near the siren and to make experiments on the focusing of intense 


sound waves in air. 



APPENDIX A 


AN AHECHOIC CHAMBER FOR USE AT HIGH AUDIO FREQUENCIES 
A. Introduction 

Careful freefield measurements of the performance of the horn 
drivers required the construction of an emechoic chamber. The requirement 
that the chamber be anecholc only for hi^ frequencies, say, frequencies 
above approximately 1 kHz, considerably simplified the acoustical treat- 
ment needed. The chamber vas to be constructed at absolute miniimnn cost, 
if possible utilizing Navy shipboard sovind insulation materials available 
as government surplus. While many construction details were forced upon 
us by budget limitations, the results may be of possible use to others. 

We found the literatvire on anecholc chambers to be surprisingly limited. 

A room on the second floor of Bxiildlng lA at Balcones Research 
Center, The Tfalversity of Texas at Austin, was available for conversion to 
an anechoio chamber. The inside dimensions were 5.0 m (hei^t) by 5*5 ® 
(width) by 7.6 m (length). The walls were of glazed brick and the ceiling 
was made of prestressed concrete U-channels. (The inverted U shape gave 
the appearance of a beamed celling. ) The floor was pointed concrete and 
sloped at a grade. Thus the hei^t was actually 5*0 ™ ®t the door and 
5.5 m at the end wall. 

We first sou^^t cm estimate of the sotind absorbing capabilities 
of the materials available to us. Althou^ financial considerations 
ultimately forced the use of the shipboard Insulation, some additional 
Justification for our selection was desired. To explain the tests per- 
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formed, we first define some terms commonly encoimtered in room acoustics. 
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Sabine absorption coefficient iiiis coefficient is the 

ratio of the sound energy absorbed by a naterial to the incident sound 
energy. For a given aaterial, the ratio is a function of frequency Euid 
the angle of incidence. The coefficient is usually measured by o<mtparing 
the reverberation times for a reverberation chamber with and vithout 
samples of the absorbing material present. 

Reverberation time T . ISie reverberation time is the time 
arequired for a sound to decay by 6o dB. This time is a function of 
and therefore is also a function of frequency. For a ro<mi that is not sc 
large that atmospheric absorption becomes important, T is given in secc'-is 
by 


^ _ (0.l6l)V 

s 57“ 


(A-1) 


Sab 


where V is the room volume, S is the room's surface area, and is the 

average value of for the room surfaces. All variables in Eq. A-1 are 

to be expressed in metric units. 

Statistical (energy) soimd absorption coefficient a^ . This 

coefficient is the ratio of absorbed-to- incident sound energy for a 
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perfectly diffuse incident sound field. Morse and Ingard state that the 

mAvinnun value Of a theoretically possible for a flat surface is about 
s 

0.96. By contrast, the value of a for a wedge-covered anechoic chamber 
wall may equal 0.996. The cutoff frequency of an anechoic chamber is 
defined as the frequency at which a drops to 0.99 (i»«*^ ® sovind reduc- 
tion of 20 dB for one reflection). Values of obtained by reverbera- 
tion room measurements usually exceed the true value of a ; sometimes the 
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excess Is 20 to 50^. For highly absorbing surfaces, aeasured values of 
°*Sab exceed 1 (Impossible by concervation of energy) by 20 to 50^. 

The excesses are attributed to difftractiou around the sample and diffi- 
culties in producing a truly diffuse field. 

Reverberation measurements vere made in a similar room using a 
"slap board" (an enlarged version of the "hcxrse vhip" device used in 
orchestras) as a sound source, a B8dC l/2 in. condenser microphone, and the 
Hevlett-Faclcard 558 OA analyzer (see chapter V) functioning as a narrowband 
graphic level recorder. The reverberation time in a 100 Hz bandwidth was 
recorded for the empty room for various frequencies. Several xypes of 
absorptive materials vere sequentially added to the chamber and the 
reverberation measurements vere repeated. The absorptive materials 
Included Owens Corning types R-11 and R-19 industrial fiberglass insula- 
tion, Oweus Coming type No. 705 acoustic insulation, and the shipboard 
fiberglass. The latter consisted of a 1') cm thick fiberglass mat with a 
perforated plastic masking on one side. The mats were 6l x 91 cm in size. 
The results of the reveroeratlon measurements overestimated by as 

much as 500^* That is, 5^^ values approaching 5 were obtained at certain 
frequencies. We made no attempt to diffuse the sound field, as is done in 
actual reverberation rooms. The presence of room modes plus the diffrac- 
tion grating effect of the besuned ceiling may have caused the poor 
results. 

Reflection coefflcien'*' measurements (normal Incidence) were made 
in the progressive wave tube described in chapter V. The equipment used 
was that shown in Pig. 5“1> with the exception that the cmechoic wedge 
termination was replaced by a sample of the test material backed by a 
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solid sluminun plxig< The results for tvo layers (10 cm) of shlpbocu'd 
fiberglass, which was the thickness eventually used for the chaaber, were 
as follcws; 


TABI£ A«1 

RfinSCnON TESTS FOR TWO lAYERS OF SHII^OARD FIBERQIA8S 
(PlASnC MASKING TCWABD WALL) 


Pfequenvy (Ht) 

200 

400 

800 

Ik 

2k 

3k 

5k 

Pressure 

Reflection 

Coefficient 

0.63 

0.35 

0.19 

0.'^3 

0.17 

0.05 

0.1 

Normal 

Incidence 

(Energy) 

Absorption 

Coefficient 

0.60 

0.89 

0.96 

0.95 

0.97 

1.0 

0.99 


Impedance tube measurements were not performed, although they might have 
improved our estimates gsreatly. In any case, for the ^quency range of 
interest, we rated the shipboard fiberglass as the best absorber among the 
materials considered. 

6. Construction 

A small section at the front of the room (J.O m (hel^t) by 
5.5 m (width) by 12.0 m (length)) was converted to an alcove to contain 
the required electronic equipment and personnel. Two separate walls of 
1/2 in. gypsum board nailed to (separate) 2 in. x ^ in. studs were 
installed between the anechoic chamber proper and che equipment room. 'Rie 
space between the walls was filled with fiberglass. The partition provided 
sufficient acoustic isolation for the personnel at high frequencies. 
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altVioxtgh nu dt>taUov t-oxu«d tranwalealon data wre> takfu. Km t\v tx- Utr 
ohfUBber tYo» th© equipaent alovve was through tvo polld wx'rxtrii iU'ura, unr 
f'or each wail« The doors h«d weather- at j ipped aeaXe. K ir^'tr tc-d ovinu c- 
tiona to equlproent uaed In the chamber were made via an erfert "tVed- 
through panel." 

The fiberglass panels ''ere attached to the walla and celling by 
using Tactoo aelf-adhe fine Insulation hangers made by AOM Industries. 

These hangers were essentially nails* 11 cm long* with self-adhest/e base 
plates. The insulation was impaled on the nails. Hroall frictlcx-flt 
retainers held the insulation tightly In place. Use of the hMigers 
simplified ccnstructlon greatly. Several alternative mounting methods 
were considered, but all would have required more installation time. 

The poured concrete floor was treated as follows. Kloor drains 
we^'e plugged and caulked to prevent sound transmission to other i^rts of 
the building. Water condensing on the floor would eventually Ixave soaked 
any fiberglass In direct contact j therefore heavy plastic sheeting was cut 
down, followed by a layer of soft rubber, 2 to 5 cm thick. 'Ibe rubber, 
which was also government surplus material, would serve to 'chanlcall;. 
isolate vibrating sources from the building. Finally, two layers v>f the 
fiberglass p*xnels were simply laid on the rubber. The panels were 
arranged to leave opt i a 0.6 m wide walkway along the room diagonal. The 
walkway was covered while experiments were under way by « roll of (wens 
Corning R-19 fiberglass. 

An overhead monorail support for microphones was suspended from 
the celJliig along the room diagonal. The total unobstructed path length 
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was 9*1 n>* A window-type air conditioner helped maintain a constant 
comfortable temperature emd humidity. 

C. Performance Test 

The traditional performance test for anecholc chambers is the 
Inverse square law test. The following test was performed using a bare 
horn driver (no horn), which acted as an approximately hemispherical 
source, except at the hlc^est frequencies used. The propagation path was 
the room diagonal, and the source was placed approximately 1.^ m above the 
floor. The microphone was carried by the monorail system. The measure- 
ments were accurate to ±0.5 dB. The propagation data are plotted in 
Fig. A and compared against the theoretical predictions for spherical 
spreading. Air absorption was negligible over such short distances and 
was not included in the calculations. We note that deviations froip 
spherical spreading were within ±3.0 dB for propagation dlsteuices out to 
5 m for frequencies between 25O Hz end 16 kHz. At the final propsigation 
(6.^ m), the microphone was pushed against the chamber corner. The effect 
of reflections near the corner sure evl^’ent in the data. The results are a 
considerable improvement over those reported by Bedell (see, for example. 
Ref. C) who also used flat absorbers rather than wedges, which are more 
effective but are also more expensive. Bedell noted deviations of ±5 
out to 1.5 m and ±5 dB out to 5 ® for bhe frequency range 500 Hz to 5 hHz. 

Our check was performed under steady state conditions. 
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Hedegaard suggests a method by which anccholc chamber perfor- 
mance may be improved. Steuidlng waves in a chamber may be reduced to 
negligible amplitudes by driving the device under test with a weurble-tone 
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and analyzing the signal with a slave filter. The changing frequency 
components are variahly phase shifted on reflection at the wall and 
therefore cemcel and ccxnpensate each other. In other voris, a steady 
state condition is never attained. The swept frequency resijonse tests of 
the drivers In the chamber were influenced by this effect, if we assume 
the sweep speed was fast enou^ to hinder mode formation. 

D. Chamber Instability 

The insulation hangers on the chamber ceiling bec^ loosening 
at many points during the late spring and svimmer of 1976. The driver 
tests had already been c<mipleted by this time. Consxxltations with the 
hanger manufacturer and distributor pxroduced several possible reasons for 
the problem but the exact cause was not pin-pointed. Because the hangers 
were designed for heating and cooling ducts used at temperatures up to 
180 * P, thermal problems frcm the aim-heated roof were not likely. It was 
concluded that either chemicals leaching from the pre stressed concrete, 
Insufficient cleaning of the ceiling surface, or a surface too porous for 
good adhesion was to blame. Since the desired driver measurements had 
already been made by this time, no serious hindrance to the i^ogress of 
the outdoor experiments was ceused. At any rate, we took consolation in 
some words expressed by Alflred E. Sommer: "Ihe best equipment is one 

that falls apart after the last experiment; any effort to make it more 
durable is a waste of time and effort." 


AFPBIDIX B 


RADIATION OF HIGH INTBISm SOUND WAVES FROM A HORN 

In this appendix ve discuss the effects of diffraction and 
nonlinear propagation distortion on the radiation of high Intensity sound 
from a circular horn. Our Interest In the subject vas triggered by two 
stimuli. First, vhlle Investigating the performance of various horn 
driver and horn combinations, we were surprised to observe high amplitude 
waveforms that did not seem to distort according to nonlinear theory 
predictions. Second, Ostrovsky and Sutin have recently published a 
series of papers dealing with the effects of diffraction on nonlinear 
propagation We briefly cite the conclusions of lineaur theory 
for diffraction by a circular horn, describe the results of a diffraction 
experiment performed at high sound levels, and finally note the predic- 
tions which follow for the experiment from the theory of Ostrovsky 
and Sutin. 

It was stated eaurller in this thesis that a ccxnmon model for 
radiation from a horn is the assvimption that a "slug" of air oscillates 
uniformly at the horn mouth. This model is equivalent to the model of 
a vibrating circular piston mounted in a plane infinite baffle, probably 
because the back wave is suppressed in both cases. The circular piston 
model is in turn mathematically equivalent to the diffraction of a plane 
wave by a circulsu: aperture in a plane screen. The solution for the 
diffraction of a plane wave may be described by a Green's function 
formulation as (see, for example. Ref. 9 ) 
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P = 



ii(t-r'/c^) 

r* 


dS 


(B-1) 


where p Is the acoustic pressure, u Is the piston velocity or the 
particle velocity wave incident upon the aperture, t is time, the dot 
indicates a time derivative, and the Integral is over the aperture 
area. Here r' is the distance from an infinitesimal area segment in 
the aperture to the axial field point. The exeu:t solution for the 
pressure at the axial position x is 

f ' Vo ( "(* - r) - - f;)] ' 


n — 2 

where r= Vx +a and a is the aperture radius. A result of Eq,. B-2 is 

that the axial pressure signal in the farfleld (Fraunhofer zone) is 
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the time derivative of the source wave. The point at which farfield 

behavior begins is known as the Rayleigh distance where 

is the piston or aperture eu:ea and X. is the wavelength] . 

A more realistic description of Taora radiation may be based 

on the Helmholtz-Kirchhoff diffraction theory for spherical waves. The 

waves radiated by a born actucLUy have a finite curvature. Thus an 

incident spherical wave is diffracted by a plane circular aperture, as 

shown in Fig. B-1. Here r is the radius of curvature of the waveform, 

o 

Q is the observation point, I is the wavefront Just filling the aperture, 
and and are the distances shown. The Kirchhoff diffraction integral 
(see, for example. Ref. 14) may be evaluated for this situation to yield 



1 + cos 9 



p(Q,t) 


• (B-3) 



FIGURE B-1 

DIFFRACTION OF SPHERICAL WAVE 
BY AN APERTURE 




The solution is interpreted as follows: The farfield pressure signed 

is the sum of a direct wave traveling edong the axis emd a phase inverted 
diffracted wave frcan the aperture rim. The diffracted wave is scaled by 
the obliquity factor (l+cos cp)/2. This result is only an apparent one 
formed by the addition of waves frcm all points in the aperture. In 
the farfield the obliqxiity factor is approximately equal to one. The 
farfield axial pressure then reduces to 


■farfield 


r + r. 
o 1 



^ p(t) 


(B-4) 


where p(t) is the incident wave in the aperture and T=t-(r^/c^) is the 
retarded time. Thus the earlier qualitative result for piston radiation, 
the wave differentiation in the feurfield, is carried over. We have 
ignored nonlinearity and atmospheric absorption in the preceding 
discussion. The omission of absorption is reasonable for audio 
frequencies over the short distances cohsldered here. The role of 
nonlineeurity is dlscusued presently. 

The experimental arrangement for the measurement of the 
waveforms on the axis of a circular horn is shown in Fig. B-2. (The 
photograph was taken in the anechoic chamber at Cockrell Hall, The 
University of Texas at Austin; the particular set of measurements 
described here were actually made outdoors). The horn vas the unit 
described in chapter V. The sovind source was a JBL 2k?0 horn driver 
(2.^ cm throat diameter). A flared connector Joined the horn and 
driver. The waveforms received are shown in Fig. The SPL on axis 

at the horn mouth was 132 dB. The waves within the hom were sufficiently 
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FIGURE B-2 

APPARATUS TO MEASURE ACOUSTIC SIGNAL 
ON THE AXIS OF A HORN 


DISTANCE OSCILLOSCOPE 
(m) GAIN (dB) 




0 



FIGURE B-3 

OSCILLOGRAMS OF RECEIVED ACOUSTIC SIGNAL 
versus AXIAL DISTANCE FOR f = 7.36 kHz 


Previously roportod by Thoobold, Wobstor, 
ond Blacksteek (1976). Rof. 48 
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intense to have suffered substantial nonlinear distortion by the time 
they reached the horn mouth. Further distortion would have tended to 
sharpen the wavefoias into a sawtooth shape. Instead, the propagating 
waves exhibited an asymmetrical shape. The waveform at a distance of 
3 m strongly resembled the time derivative of the mouth waveform 
(distance=0) . Unfortunately, the mouth waveform also contained a 
contribution from the diffracted wave such that the soiirce wave could 
not, in general, be observed directly. 

Because the boundary conditions for the Kirchhoff solution 
(Eq,. B-3) are not exact, i.e., the solution cannot reproduce the 
assumed conditions in the aperture plane, then the solution can only be 
approximately correct near the aperture. We may get an indication of 
the expected behavior in the aperture for large source wave curvatures, 
r^, by examining the plane wave solution (Eq. B-2). The exact plane 
wave solution on- axis in the aperture is 

- f-) • 

o 

We note that the horn radius for our experiment was 7*8 c® afiti "the time 
delay corresponding to a/c^ was very nearly (3/2) T, where T is the 
period. By periodicity, we expect that the pressure amplitude on-axis 
at the mouth was twice the amplitude of the incident wave. The wave 
shape is expected to remain unchanged. We assume that the source wave 
calculated here may be used for the spherical diffraction problem. 

As a rough quemtitatlve check on the applicability of the above 


euialysis, we calculated t..e peak- to- peak amplitude of the source wave 
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based on the experimental farfleld signal. This peurtlcular course vas 
followed because it is much simpler to numerically integrate than to 
differentiate. The relation anployed is 

r + rj^ c - 

^aperture ” r r^ - r, J ^farfield (®~^) 


The farfleld signal at a distance of 3*0 m was integrated and scaled as 

in Eq,. B-6 to predict a mouth waveform of a peak- to- peak pressure ampli- 
2 

tude 170 N/m . The pressure amplitude measured at the horn mouth was 
2 

240 N/m (p-p). We have shown that the actual incident wave should 
have an amplitude of about half that of the wave at the mouth, or 

O 

120 N/m*" (p-p). The agreement is within 3 dB. While the error is 
fedrly l€u:ge, it is not too surprising, considering the variablility 
in the horn mouth pressure distribution noted in chapter V, and that 
the microphone was aligned with the horn by eye and the measurements 
were made outdoors. Nonlinearity should not have played a Icurge role 
following diffraction because the amplitude was sharply reduced from 
an orlginadly marginal level by spherical spreading. 

We now briefly discuss a few points from Ostrovsky's work on 
diffraction problems in nonlinear acoustics. A geometry shown below 
is assumed. 




FIGURE B-4 


DIFFRACTION FROM A HORN (after Ostrovsky and Sutin^^) 
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Consider a horn radiating finite amplitude sound generated by a vibrating 

piston in the horn throat. The source wave has the curvature R . The 

s 

distortion variable, expressed in stretched coordinates, for propagation 
through the horn to the mouth is (see, for example. Ref. 4) 

where A is the throat area A(x) is the cross sectional area at the 
o 

position X, eind quasipleme wave propagation through the horn is assumed. 
Ostrovsky and Sutin base xheir analysis on the hypothesis that, for 
certain geometries, frequencies, and source levels, the diffraction 
and nonlinear distortion stages may be separated. The stated condition 
for separability is 


dx 


iB-7) 


Pe « > vB-8) 

where a is defined in Eq,. B-7 and e is measured at the horn throat. 

(Shock formation is precluded here. For cases involving the diffraction 
of shock waves, see Ref. 74.) We offer a physical reason in support of 
the idea of separability. The only nonlinear effect considered here, 
namely, propagation distortion, is an effect cumulative with distance. 
Diffraction, on the other hand, is a localized phenomena. If the wave 
is not too strong, the propagation distortion in the Fresnel zone will 
be negligible and the diffraction and distortion may be considered 
separately. Van Buren and Breazeale' have utilized a parsdlel. hypothesis 



in describing the reflection of short fi, Ite amplitude sound pulses 
from rigid boundeuries. 

The values of the parameters for the measured waveforms shown 

in Fig. B-3 relevant to the theory of Ostrovsky and Sutin were as 

fellows: 3e-i2.&<10 ^ , o=0.5 (calculated value), R =23 cm, and 

8 

R =4l cm. Under these conditions, the diffraction and nonlinear 
o 

distortion stages are predicted to be sepsu^able. The solution given by 
Ostrovsky and Sutin reduces to Eq, B-5, where the obliquity factor is 
equal to one. A second solution applicable for observation points off 
the horn axis is also given. We have not experimentally investigated 
this case for lack of time. 
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